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The  relationship  between  the  product  energy  distribution  of 
a  bimolecular  reaction  and  the  energy-dependent  cross  section  of  the 
reverse  reaction  is  examined.   In  the  case  where  the  forward  reaction 
is  associative  ionization 

A  +  B  =  AB  +  e 
the  reverse  reaction  will  be  dissociative  recombination 

AB'^  +  e~  =  A  +  B  . 
It  is  pointed  out  that  experimental  measurements  of  the  rate  coeffi- 
cient for  dissociative  recombination  may  be  in  error  due  to  electron 
heating  during  the  recombination  process  or  due  to  a  source  of  non- 
thermal electrons.   The  conditions  under  which  this  may  be  significant 
are  discussed.   An  expression  for  the  temperature  dependence  of  the 
dissociative  recombination  rate  coefficient  is  derived  using  the  prin- 
ciple of  microscopic  reversibility  and  an  expression  for  tlie  associa- 
tive ionization  cross  section  which  is  derived  from  the  classical 
solution  to  the  atom-atom  ionization  cross  section.   The  result  fits 

iv 


the  experimental  results  and  is  obtained  without  making  any  detailed 
assumptions  about  the  nature  of  the  A  -  B  repulsive  potential  curves. 

A  general  theory  for  product  energy  distributions  is  developed 
by  modifying  the  statistical  phase  space  theory.   The  assumption  is 
introduced  that  the  probability  of  product  formation  depends  on  the 
momentum  ti-ansferred  in  the  reactive  collision.   Minimum  monentura  trans- 
fer in  the  center-of-mass  coordinate  system  is  assumed  to  gi^'e  the 
greatest  probabilitj'  for  reactive  scattering.   The  effect  of  the  momen- 
tum transfer  modification  is  shown  by  calculations  of  the  reaction  cross 
section  and  probability  of  reaction  as  a  function  of  impact  parameter 
for  the  reaction  K  +  HCl  =  KCl  +  H.   These  calculations  are  compared 
with  results  calculated  from  the  unmodified  phase  space  theorj-  and  with 
publislied  molecular  beam  results.   The  raodified  phase  space  theorj' 
gives  results  closer  to  the  experimental  values  tlian  the  unnodified 
theory.   The  general  implications  of  the  theory  are  discussed. 
Particula2-ly  interesting  is  the  possibilitj'  that  such  a  modified 
theory  may  lead  to  more  realistic  vibrational  energy  distributions 
than  the  pi'esent  theory. 


I.   INTRODUCTION 

The  product  energy  distribution  for  bimolecular  reactions  of 
the  type 

A  +  B  =  C  +  D  (1) 

has  been  studied  by  a  number  of  workers  for  a  number  of  specific  cases. 
The  theoretical  problem  is  generally  too  difficult  to  solve  using  exact 
quantum  mechanical  methods,  so  a  number  of  semiclassical ,  classical, 
and  statistical  approximations  have  been  used  to  obtain  numerical  solu- 
tions. 

Most  of  the  experimental  work  and  theoretical  attention  has  been 
directed  to  exothermic  reactions  [1-6] .   Exothermic  reactions  often 
have  large  cross  sections,  a  large  amount  of  energy  to  distribute  among 
the  product  states,  and  frequently  yield  energetically  excited  species 
which  are  readily  detectable  by  spectroscopic  or  chemical  means. 

However,  endothermic  reactions  have  useful  features  which  can 
be  exploited  too.   The  small  number  of  energetically  accessible  product 
states  near  the  threshold  energy  makes  it  possible  to  describe  the 
product  states  more  completely  than  in  the  exothermic  case.   Also,  the 
use  of  the  principle  of  microscopic  reversibility  makes  it  possible  to 
apply  the  results  of  endothermic  reaction  measurements  and  calculations 
to  the  study  of  the  corresponding  exothermic  reaction  [7]. 


The  experimental  determination  of  product  energy  distributions 
is  limited  by  the  thermal  distribution  of  reactant  states  and  by  the 
rapid  thermal ization  of  product  states.   This  limitation  may  be  over- 
come by  employing  molecular  beams  with  velocity  selectors  for  both 
reactants  and  products.   Also,  the  use  of  coincidence  techniques  enables 
experimentalists  to  measure  product  energy  distributions  directly. 
Measurements  have  been  limited  by  the  need  to  obtain  a  sufficiently 
large  signal.   Recently,  Zare  and  co-workers  [8]  have  pioneered  the  use 
of  laser-induced  fluorescence  to  measure  the  product  energy  distribu- 
tion in  molecular  beams.   Other  experimental  determinations  of  product 
energy  distribution  in  molecular  beams  have  been  reviewed  by  Kinsey  [5]. 

Another  method  of  determining  product  energy  distributions  is 
the  measurement  of  chemiluminescence.   The  use  of  this  technique  has 
been  reviewed  recently  by  Carrington  and  Garvin  [2],  Polanyi  [9],  and 
Carrington  and  Polanyi  [4].   Anlauf  et  al .  [10]  have  shown  how  the 
experimental  data  of  a  chemiluminescent  exothermic  reaction  may  be  used 
to  calculate  the  rate  of  the  reverse  endothermic  reaction. 

Each  of  these  techniques  has  limited  applications  to  endothermic 
reactions.   The  small  cross  section  of  most  endothermic  reactions  limits 
their  observation  in  molecular  beams,  and  the  limited  amount  of  energy 
available  rarely  results  in  chemiluminescence. 

Another  method  that  may  be  used  in  certain  cases  is  the  measure- 
ment of  product  temperatures  by  spectroscopic  means.   This  requires  that 
the  product  states  thermalize  slowly  with  respect  to  the  time  of  the 
measurement. 


This  condition  is  met  by  the  associative  ionization  jjroccss 


A  +  B  =  Ab"*"  +  e"  (2) 


which  has  been  observed  in  plasmas  [11-19] ,  shock-tube  investigations 
[20-23]  and  molecular  beams  [24] .   This  reaction  is  endothermic  unless 
one  of  the  reactants  is  electronically  excited,  in  which  case  it  may  be 
considered  as  a  special  type  of  Penning  ionization  [25] .   The  reverse 
reaction,  known  as  dissociative  recombination, 


Ab"^  +  e   =  a  +  B  (3) 


is  very  important  in  plasma  recombination  studies,  due  to  this  large 
cross  section,  and  has  been  studied  extensively  [26-50]. 
The  associative  ionization  reaction 

N  +  0  =  NO"^  +  e~  (4) 

is  endothermic  by  270.  kj  mol   .   It  is  important  as  a  potential  source 
of  the  NO   ion  because  of  its  low  energy  threshold  and  has  been  studied 
in  shock  tubes  [20-22]  and  in  connection  with  ballistic  reentry  [51,52]. 
Due  to  the  light  mass  of  the  electron,  it  will  take  a  relatively  long 
time  for  the  product  electrons  to  reach  the  temperature  of  the  surround- 
ing gas  if  they  are  produced  at  a  much  higher  or  lower  temperature  than 
the  ions.   Thus,  it  sliould  be  possible,  by  measuring  the  temperature  of 
the  electrons  produced,  to  determine  the  product  energy  distribution  of 
this  endothermic  reaction. 

It  is  necessary  to  calculate  the  product  energy  distribution 
for  this  reaction  in  order  to  obtain  accurate  experimental  rate  measure- 
ments of  processes  involving  electrons  in  systems  where  associative 


ionization  is  an  important  source  of  electrons.   The  mass  disparity 

between  electrons  on  the  one  hand  and  positive  ions  and  neutral  species 

on  the  other  results  in  slow  momentum  transfer  between  the  electrons 

and  the  heavier  particles,  although  particles  of  the  same  mass  will 

transfer  momentum  efficiently.   Thus,  if  the  electrons  are  formed  with 

more  or  less  kinetic  energy  than  the  average  kinetic  energy  of  the  ions 

and  neutral  species,  a  two-temperature  plasma  will  be  formed,  with  the 

electrons  having  an  energy  distribution  corresponding  to  an  electron 

temperature  T   and  the  more  massive  particles  having  an  energy  distri- 
e 

bution  corresponding  to  another  temperature  T  .   This  may  affect  exper- 
imental measurements  of  processes  involving  electrons  because  many  such 
processes  have  a  strong  dependence  on  electron  energy  and  because  many 
of  the  procedures  for  determining  electron  density  (the  Langmuir  probe 
and  microwave  diagnostics)  are  dependent  on  electron  temperature. 

We  have  reported  previous  calculations  [53]  of  the  product 
energy  distribution  for  associative  ionization.   In  these  calculations, 
it  was  assumed  that  the  colliding  atoms  were  hard  spheres  and  that 

reaction  took  place  whenever  the  line-of-centers  energy  E   exceeded 

^•'   Ic 

the  threshold  energy  E   for  the  reaction, 

E.,   =  E.  (1-b^/R^)  <■   E*  (5) 

Ic     1 

where  E.  is  the  initial  relative  kinetic  energy,  b  is  the  initial 
impact  parameter,  and  R  is  the  hard-sphere  collision  radius.   It  was 
pointed  out  that  almost  all  of  the  relative  orbital  angular  momentum 
of  the  reactants  was  transformed  into  i-otational  angulai-  mom.entum  of 
the  product  molecular  ion.   Also,  it  was  assumed  that  the  probability 


of  the  production  of  product  in  vibrational  level  v=  i,  given  a  total 
energy  available  for  distribution  between  vibrational  energy  and  elec- 
tron kinetic  energy,  was  independent  of  vibrational  state  and  directly 
related  to  the  probability  of  producing  electrons  with  energy 


e     =  €   -  ihv, 
e    n 


P(v=i)    g(en-"^^> 


(6) 


P(v=0)      g(G^) 


It  was  assumed  that  g{£   )  was  a  function  of  the  form 


g(€  )  =  (€  )  ""  (7) 

e      e 


where  n  could  be  established  to  be  between  0.5  and  1.5  on  the  basis  of 
experimental  measurements  of  dissociative  recombination  rate  coeffi- 
cients. 

The  results  of  these  calculations  showed  that  electrons  v.ere 
produced  at  a  lower  temperature  than  those  of  the  neutral  species. 
This  suggested  that  care  should  be  taken  when  measuring  dissociative 
recombination  rates  of  ions  produced  by  associative  ionization.   If 
the  electrons  and  neutral  species  have  not  yet  thermalized,  the 
results  can  be  erroneous  if  the  analytical  technique  used  is  sensitive 
to  electron  temperature. 

This  calculation  also  raises  the  question  of  the  electron 
thermalization  rates  in  a  recombining  plasma.   If  dissociative  recombi- 
nation is  more  rapid  for  lower  energy  electrons,  higher  energy  electrons 
will  remain  in  a  recombining  plasma  resulting  in  an  increase  in  the 
electron  temperatui-e.   The  possible  magnitude  of  this  effect  is  out- 
lined in  tlie  laext  section. 


II.    ELECTRON  THERJ.IALIZATION  IN 
A  RECOMBINATION-DOMINATED  PLASMA 


A.   Basic  Equations 


In  order  to  understand  the  effect  of  the  recombination  on  the 
thermalization  time  of  the  electrons  and  on  the  measured  recombination 
rate  coefficient,  it  is  necessary  to  construct  a  model  of  electron 
temperature  changes  in  a  recombination-dominated  plasma.   First,  it  is 
assumed  that  diffusion  and  radiative  recombination  are  negligible  com- 
pared to  dissociative  recombination  and  that  there  is  no  source  of 
electrons  and  ions  in  the  recombining  plasma.   Then  the  rate  of  electron 
loss  may  be  exjaresscd  as 


dn 


e  (8) 


=  -  a  n  n 
dt       R  e  + 


where  n   is  the  electron  density,  n     is  the  dissociative  recombination 
e  R 

rate  coefficient,  n   is  the  concentration  of  positive  ions,  and  t  is 
the  time.   If  it  is  assumed  that  n  =  n  ,  then 


dn         2 

=  -  cv  n  (9) 

dt       R  e 


and,  if  g'   is  constant,  then  at  time  t, 
R 

o 

„£  ^  n°c^t   +   1  (10) 

n      e  R 
e 


o 

where  n   is  the  electron  concentration  at  t = 0.   Now,  if  the  relaxation 
e 

time  T  is  defined  as 

1    o 

-  =  n^cv^  (11) 

and  the  reduced  electron  density  N  as 

N  =  ^  ,  (12) 

n 
e 


then  Equation  (12)  can  be  expressed  as 


N   =  1  +  t/T  .  (13) 


So  far,  we  have  treated  the  dissociative  recombination  rate 

coefficient  as  if  it  were  a  constant.   In  fact,  it  is  a  function  of 

the  electron  temperature  and,  like  the  electron  temperature,  •.-.ill  vary 

with  time.   To  consider  this  more  general  case,  let  a     be  the  recombi- 

R 

nation  rate  coefficient  when  T  =T,  and  let 

e   A 


\^e 

A  dimensionless  differential  equation  can  be  obtained  bj' 

o  2  A 
dividing  both  sides  of  Equation  (9)  by  (n  )  Op  to  obtain 

d(n  /'n°)      ,a.      n  ,2  (15) 

e   e      f  R\   I    e> 


or 


l^j  {-o) 


d(t/Tj       - 

R    e 


dN      ^2 

de=   -'"-^ 


.  ,  A  ,  o 

where   8  =  t/T  a  -  &„/a„ ,   and  N  =  n  / n  . 

A       R   R  e   e 


The   total    electron   energy   per  volume,    E^^_^,    will    be 


(17) 


tot  e   e 

where  6   is  the  average  energy  of  the  electrons.   Then 
e 

dE    .  =  (-J2l)^,   +  f-^^de_  (18) 


'tot    V  Sn^  y   e    X  g^ 


~e 


or 


Then, 


dE     =  f  dn  +  n  df 
tot   ^e   e    e  e 


dE,  .      dn       d^ 


(19) 


dt   "  "^  dt     e  dt 


Rearranging,  and  substituting  Equation  (9), 


"!£  =  i  f^  .  c  n^?  )  .  <21) 


dt    n   V  dt      R  e  e 
e 


The  rate  of  change  of  the  total  electron  energy,  dE^^^'dt, 
will  depend  on  the  electron  loss  and  electron  production  processes 
which  are  important,  and  on  energy  exchange  between  electrons  and 
other  species  in  the  system.   If  we  assume  that  electrons  are  not 
being  produced  and  that  the  dominant  electron  loss  process  is  disso- 
ciative recombination,  then  only  three  terms  will  contribute  to  the 
total  electron  energy:   dissociative  recombination,  elastic  collisions, 
and  inelastic  collisions.   It  has  been  shown  that  the  contribution 
from  inelastic  collisions  may  be  neglected  for  electrons  with  energies 
less  than  the  vibrational  quanta  of  the  neutral  gas  molecules  [54,55]. 
The  total  energy  change  is  then  given  by 

dE^  ^    dE  ^    dE 

tot  _    el     RC  (22) 

dt   "   dt   '*'   dt   ■ 


Petschek  and  Cyron  [56]  considered  energy  transfer  between 
electrons  and  neutral  particles  and  obtained  the  following  expression 
for  the  change  in  electron  energy  due  to  elastic  collisions,  dE^^/dt, 


^  ,T.     .  ,,    ^3/2 


""el  _  4  %  _  ^   C_A_   ,V'4  -.  ^  (23) 


^F-=  -^  VA^rr  \T--   V^3  -e 

TT-   A  e 

where  m   is  the  electron  mass,  m^  is  the  mass  of  the  neutral  particles, 
e  ^ 

n   is  the  neutral  gas  density,  and  o^^^   is  the  momentum  transfer  cross 
section  for  electron-neutral  collisions.   The  momentum  transfer  cross 
section  has  been  measured  for  a  number  of  systems  [57-66]  and  is  gen- 
erally a  function  of  the  electron  energy.   However,  here  it  is  assumed 
that  the  momentum  transfer  cross  section  is  a  constant.   This  approx- 
imation is  valid  only  if  a^^^   changes  very  slowly  as  a  function  of 
energy,  or  if  the  electron  energy  changes  slowly.   Experimental  mea- 
surements and  energy  dependences  of  the  momentum  transfer  cross  sec- 
tions of  electrons  in  several  gases  are  discussed  in  Appendix  A.   It  is 
clear  that  no  simple  function  describes  the  energy  dependence  of  all  the 
cross  sections  discussed.   It  is  apparent,  however,  that  for  molecular 
gases,  such  as  0^  and  N^ ,  that  the  momentum  transfer  cross  section  is 
fairly  constant  at  low  energies. 

It  is  a  little  harder  to  justify  this  approximation  with  rare 
gases,  on  account  of  the  pronounced  minimum  which  appears  in  their 
momentum  transfer  cross  section.   The  calculations  of  recombination 
heating  to  be  discussed  later  show,  however,  that  heating  effects  are 
likely  to  be  important  only  when  the  electron  temperature  is  fairly 
close  to  the  gas  ter.perature.   Thus  the  approximation  of  letting  the 


10 

momentum  transfer  cross  section  be  constant  can  be  justified,  providing 

that  a  value  is  always  chosen  that  is  consistent  with  the  neutral  gas 

temperature.   Notice  that  tlie  total  electron  energy  will  either  decrease 

or  increase  due  to  elastic  collisions,  depending  on  whether  the  electron 

temperature  is  higher  or  lower  than  the  neutral  temperature.   Wlien 

T  =  T, ,  then  dE  , /dt  =  0. 
e   A         el 

The  total  electron  energy  will  also  change  as  the  result  of 
dissociative  recombination.   The  electron  energy  change  due  to  recom- 
bination, dE  /dt,  is  given  by 
RC 

dE         dn 

^^   €„.  -^  (24) 


dt     RC  dt 


or 


— ^  ^     -a  f     n^  (25) 

dt       R-RC  e 

where  the  average  energy  of  the  recombining  electrons,  (^      ,    is  given  by 

RC 

03 

J     6aj,(6)€  exp  (-6/kT^)de 

'°    ! (26) 


^RC     «» 

f   CT^(t)€  exT  (-6/kT  )df 
^         R  e 

o 

where  a   (c)  is  the  energy-dependent  cross  section  for  recombination. 
R 

If  a  (O    is  of  the  form 
R 

a, (6)  =  K/6"  (27) 


then,  for  n  <  2 , 


:    T-.    (2-n)kT   .  (28) 

-RC  e 


Now,  substituting  Equations  (23)  and  (25)  into  (22)  and  then  into 
(21),  the  change  in  average  electron  energy  is  given  by 


where 


e 


3/2 


R  e 


-V 


C  =  (4M^)  im^/-n)^    n^o^^ 


Replacing  f     by  —  kT  , 
e    2    e 
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(29) 


dt  "  3k 


C(--lJi^2kT^ 
e 


3/2 


+  &n        -  kT  -  f^^ 
R  e  \2   e   RC 


(30) 


When  the  suggested  cross  section  of  Equation  (27)  is  used, 
Eqviation  (30)  becomes 


dt   3k 


C(/-lJ^2kT^ 
e 


3/2 


+  a„n  kT  (n  - -) 
R  e   e    2 


(31) 


It  can  be  seen  tliat,  for  n  =  l/2,  the  term  due  to  dissociative  recombi- 
nation vanishes,  and  the  electron  temperature  will  not  be  distui'bed 
by  the  recombination  reaction.   If,  however,  n  is  greater  than  1/2, 
the  average  energy  of  the  recombining  electrons  will  be  less  than  the 
average  energy  of  all  the  electrons.   This  will  result  in  an  increase 
in  the  average  energy  of  the  remaining  electrons  as  the  recombination 
process  proceeds. 

The  rate  coefficient,  o  ,  is  given  by 

K 


-   ,     — Y    -^   r   a    (06  exp    (-e/kT   )c]^ 

R      \m  TTkT    /     kT     >J       R^^        i     V    ^        g      - 
e        e  e    o 


Q'     = 


(32) 


The   integral    is   just   the  Laplace   transform  of   the   function   a    (6)€   so 


.        ..^      J_:f    ( J_V(2-n)(kT    )2- 
R        \m  TTkT    /      VkT    /  e 

e        e  e 


cy„  =  K 


(33) 
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or 


a     =  Kr(2-n)(— Y    (kT   )  ^"s"''^    .  (34) 

R  \m  TT/  e 


Wlien  T    =T    ,    the   recombination  rate   coefficient    is 
e        A 

,A     Kr(2-„,m*(kT)<*-">  "^' 

R  \m  TT/  A 


so  that  the  temperature  dependence  of  the  reduced  rate  coefficient 
is  given  by 

a  =  -I  =  (^)     .  (36) 

a  e 

R 

Equation  (31)  can  be  simplified  by  using  the  reduced  parameters 

N  =  n  /n° ,   9  =  t/T  and  /  =  T  /T ,  to 
e   e  *-'    e  A 

r  (2kT  )*C'        ,     c. 

^=3— A-    '"--dt  '-^l-l   Nd  (n-i)  .  (37) 

^R^e  ^R 

In  summary,  the  recombination  process  may  be  described  by  two 

coupled  differential  equations  which  are  expressed  in  terms  of  reduced 

time,  9,  electron  density,  N,  and  electron  energy,  /": 

(--n) 
|L.c^l(l^),|(n-|)r^2   "  (33^ 

wliere   the   form  of    the   solution  depends   on   tlie   constants   n    and 

^  ^     ,2m  kT.   i  n.    a 
C   .    11     -£-^Y    -i-4   .  (40) 

3      \         TT         /  o  A 

"e">A^R 
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It  should  be  noted  that  this  ti'eatment  assumes  that  there  are 

separate  Maxwell -Boltzinann  distributions  of  electron  and  neutral 

-    3  -     3 

energies,  so  that  ^  =  -  kJ      and  £     =  -  kT    This  implies  that  colli- 

sions  between  electrons  are  frequent  enough  to  maintain  the  distribu- 
tion function  as  the  average  energy  of  the  electrons,  f    ,    changes. 

e 


B.   Numerical  Examples 

In  general,  Equations  (38)  and  (39)  must  be  solved  numericallj-- 

by  an  iterative  technique.   The  apparent  reduced  rate  coefficient, 

a  =   <y  /en     may  be  taken  as  the  slope  of  a  plot  of  1/N  vs.  3  and  will  be 
R   R 

unity  for  the  "ideal"  case  where  T  =  T,  and  there  is  not  a  recombina- 
■'  e    A 

tion  heating  effect.   Effects  of  non- thermal  electrons  may  be  due  to 

electrons  being  produced  with  an  initial  T  t^   T  or  to  electrons  being 

heated  by  the  recombination  process  preferentially  removing  lower 

energy  electrons. 

A  set  of  numerical  solutions  to  Equations  (38)  and  (39)  is 

given  in  Figures  la-c  for  the  case  n=l,  C=l,  and  several  different 

values  of  the  initial  reduced  electron  temperature,/-  .   Figure  la 

indicates  how  electron  temperature  varies  with  time,  Figure  lb  shows 

the  decay  of  electron  density  as  a  function  of  time,  and  Figure  Ic  gives 

the  ratio  of  observed  to  the  ideal  thermal-to-neutral  rate  constants. 

The  curves  in  Figure  Ic  are  the  derivatives  of  the  corresponding  curves 

in  Figure  lb.   As  exiDected,  elastic  collisions  (the  first  terms  on  the 

right-hand  side  of  Equation  (33))  will  tend  to  decrease  ^when^^^  >  1 

(or  T   >  T  )  and  increase/^  v.hcn  f      <  1    (or  T  <  T.).   However,  for 
^     e    A  ^      ^o  e    A 

n  >  1/2  the  recombination  effect  (tlie  last  term  on  the  right-hand  side 
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Figure  la. 


Plot   of  Q   vs    6   foi'   different    initial    valvies   oif  :       ^^  f     ' 
b)0       =    .5,    c)(f       =    1,    d)(^       =    2,    e)  (^      =    5,    f)  /      -.10?, 

^-'ri  ^-'n  ^^n  ^^n  ^^n 
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Figure  lb. 


Plot  of  IN  vs  e  for  different  initial  values  of/  :  ^) ^  ^ 
W6^=    -5,  c)(5,=  1,  cl)(f^.  2,  e)(^^.  5,  f)t,,=    30, 
C  :.  1,  n  =  1. 
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Figure  Ic. 


Plot  of 

h)f       =    .5,  c 

C  =  1,  n  =  1. 


a  vs  G  for  different  initial  values  of  (f:   a)  (f 
""     ^   -^  ^   =  2,  e)(f   ..  5.  i)J   =  30? 


)Co  =  1.  ^')(^ 


n 

3 
2 


o 
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of  Equation  (38))  will  always  result  in  electron  heatin^'-.   For  T  <  T 

e    A 

the  heating  effects  of  recombination  and  elastic  collisions  reinforce 
each  other,  while  for  T  >  T   the  two  processes  oppose  one  another. 
Thus  "cold"  electrons  will  heat  up  faster  than  "hot"  electrons  will 
cool  down. 

The  initial  electron  temperature  will  depend  on  the  ionizing 
source  used  to  produce  the  electrons.   Very  low  electron  temperatures 
would  be  expected  from  threshold  photoionization  sources,  while  quite 
high  temperatures  could  be  obtained  from  discharge  tubes.   Obviously, 
it  would  be  easy  to  predict  the  initial  electron  temperature  in  the 
event  that  the  ionizing  source  was  photoionization  or  Penning  ionization. 

The  dimensionless  parameter,  C,  characterizes  the  relative 
rates  of  electron  thermalization  and  recombination.   Figures  2a-c  and 
3a-c  show  how  the  rate  parameters  vary  for  different  values  of  C. 
With  small  values  of  C  there  can  exist  significant  deviations  of  the 
electron  temperature  from  the  neutral  temperature.   Thus  if  electrons 
are  initially  produced  with  a  temperature  different  from  the  neutrals 
(Figures  2a-c)  or  are  heated  by  the  recombination  process  (Figures  3a-c) , 
long  term  disparities  in  the  electron  and  neutral  temperatures  can  per- 
sist (Figures  2a  and  3a) .   Because  of  tlie  temperature  dependence  of 
the  recombination  rate  coefficient,  observed  rate  constants  may  deviate 

significantly  from  those  which  would  be  measured  if  T  =  T,  (Figures  2c 

e    A 

and  3c).   However,  at  conditions  for  which  C>1,  elastic  collisions 
rapidly  tliermalize  the  electrons  to  the  neutral  temperature. 

In  general,  the  electron  temperature  and  recombination  rate 
changes  will  be  minimized  in  systems  with  a  large  momentum-transfer 
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Figure  2a. 


Plot  of(£vs 


for  different  values  of  C:   a)  C  =  . 01 , 


b)  C  =  .1,  c;  C  =  .5,  d)  C  =  1,  e)  C  =  10. Co  =  2,  n  =  1. 
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Figure  2b.   Plot  of  1/N  vs  G  for  different  values  of  C:   a)  C  =  .01, 


b)  C  r.  .1,  c)  C  =3  .5,  d)  C  =  1,  e)  C 


10, 


=  2,  n  =  1. 
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Figure  2c.   Plot  of  a  vs  8  foi^  different  values  of  C:   a)  C  =  .  01 , 


b)  C  =  .1,  c)  C  --  .5,  d)  C  =  1,  e)  C  =  10,(^ 


n  =  1. 
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Figure  3a.   Plot  of (£  vs  0  for  different  values  of  C:   a)  C  =  . 01 , 

b)  C=  .1,  c)  C=  .5,  d)  C=l,  e)  C=10,^^=l,  n=l. 
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Figure  3b.   Plot  of  1/N  vs  8  for  different  values  of  C:   a)  C 


01 


b)  C  =  .1,  c)  C  =  .5,  d)  C  =  1,  e)  C  =  10 ,  (>^   =1,  n  =  1, 
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Figure  3c.       Plot   of    a   vs    9   for   different  values   of   C:       a)    C  =    .01, 

b)    C  ^    .1,    c)    C  =    .5,    d)    C  =   1,    e)    C  =   10,  ^^  =   1,    n  =    1 
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cross  section  for  electron-neutral  collisions,  neutral  particles  witli 
small  mass,  a  small  degree  of  ionization,  and  small  recombination  rate 
coefficients.   Such  effects  are  minimized,  for  example,  by  diluting 
a  recombining  plasma  with  a  neutral  carrier  gas,  particularly  with 
a  light  noble  gas  such  as  neon  or  argon. 

As  previously  described,  if  the  energy  dependent  cross  section 
is  of  the  form 

r 
then  the  rate  constant  is  of  the  form 

1 

R    e 

Since  recombination  rate  constants  are  generally  reported  with  temper- 

-1/2      -3/2 
ature  dependences  ranging  from  T     to  T     ,  a  range  of  n  from 

1  to  2  is  suggested.   However,  foi'  the  case  n=2  integration  of  the 
cross  section  over  a  Boltzmann  distribution  leads  to  an  undefined 
result,  as  will  be  discussed  later.   If  n=l/2,  the  rate  coefficient 
is  temperature  independent  and  electron  temperature  deviations  will  not 
affect  observed  rate  constants.   Also  £     -    \        for  i\--\/2,    so  recombi- 
nation heating  will  not  occur  for  a  system  with  a  temperature  indepen- 
dent rate  constant.   But  e:-cperimental  measut'ements  of  recombination 
rate  coefficients  generally  demonstrate  an  inverse  temperature  depen- 
dence, so  n  >  1/2  is  indicated. 

The  calculated  solutions  of  Eqiiations  (38)  and  (39)  are  indi- 
cated in  Figures  4a-c  for  several  different  values  of  n.   As  expected 
from  the  previous  discussion,  larger  recopabination  heating  effects  are 
seen  for  larger  values  of  n. 


25 


2  IL 


0 


^ 


:!;«5lt»e***,i*>*wrtwifik»«»r«i»*««l»tr«W»tW»:?fc^^ 


1 


e 

Figure  4a.   Plot  of(^  vs.  6  for  different  values  of  n: 

b)  n  =  1.5,  c)  n  =:  1.0,  d)  n  =  0.5,/   =  2, 

'-'  o 


a)  n  : 
C  =  1. 


1.8, 
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Figure  4b.   Plot  of  1/N  vs  6  for  different  values  of  n:   a)  n  =  1.8, 
ta)  n  =   1.5,  c)  n  r.  1.0,  d)  n  =  0.5,  (^  .=    2,    C  =   1. 
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Figure  4c.      Plot   of    n   vs      G   for  different   values   of  n:       a)    n  =   1.8, 


b)    n 


1.5,    c)    n=    1.0,    d)    n:^   0.5,  C      =    2,    C 

o 


1. 
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C.   Otliei'  Possible  Dissociative  Recombination 
Cross  Sections 


An  appropriate  form  for  the  dissociative  recombination  cross 
section  can  be  found  by  examining  the  experimental  measurements  of  the 
temperature-dependent  recombination  rate  coefficient.   The  relationship 
between  the  two  functions  was  noted  in  Equation  (32)  above.   Measure- 
ments at  room  temperature  and  other  low  temperatures  correspond  to 

-1/2 
a  temperature  dependence  of  T     .   If  the  cross  section  is  of  the  form 

e 

o(0    ^   6   ,  then  n=l  corresponds  to  this  temperature  dependence  of  the 
rate  coefficient. 

At  higher  temperatures,  the  dissociative  recombination  rate 

-3/2 
coefficient  is  observed  to  become  proportional  to  T     ,  not  only  for 

NO   [21]  but  also  for  other  diatomic  ions  [67].   This  would  appear  to 

-2 
correspond  to  the  cross  section  a(€)  '^  £      .      However,  when  this  cross 

section  is  substituted  in  the  expression  for  the  recombination  rate 

coefficient,  Equation  (32)  and  integrated,  the  result  contains  the 

exponential  integral  with  zero  argument,  E  (0) ,  which  is  not  defined 

[68].   It  was  suggested  by  Hansen  [69]  that  the  rate  coefficient  was 

of  the  form 

a^  ^    c(kT)"    [1-exp  (-d/kT)]  (41) 

where  c  and  d  are  constants.   At  low  temperatures,  the  exponential 
term  is  negligible;  at  higher  temperatures  it  may  be  expanded 

9  ? 

d   1  /  d  \    1  /  cP 


exp  [-   _;=l-~+--  [-)    -^[~)    +...  (42) 
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and  the  rate  coefficient  is  of  the  form 

-3/2 
a^  ^   dc(kT)  (43) 

if  higher-order  terms  are  neglected. 

This  rate  coefficient  may  be  derived  from  several  different 
cross  sections.   O'f.Ialley  [70]  assumed  that  the  exponential  term 
arises  from  the  vibrational  excitation  of  the  ion  at  higher  temper- 
atures.  If  the  dissociative  recombination  cross  section  is  negligible 
for  all  vibrational  states  except  the  ground  state,  then  the  rate 
coefficient  will  depend  on  the  fraction  of  molecular  ions  in  the  zeroth 
vibrational  level.   Then  the  cross  section  for  dissociative  recombina- 
tion is 

c  =  c^e  if   V  =  0 

(44) 
Cj  =   0        if   V  =1,2,3,... 

The  cross  section  is  also  implied  in  Hansen's  derivation  of  the  rate 
coefficient.   The  rate  coefficient  corresponding  to  this  cross  section 
may  be  written  as 

cy^,  =c,(kT  )~^'^^[l-exp  (-h'^/kT.)].  (45) 

R   1   e  ■  A 

Another  possible  cross  section  is 


a  =c/a  0  ^  f  ^  w 

II    2  '^  ^ 


a       -      0  g  >  w 

II  ^ 


(4G) 


wliere  w  is  a  constant.   This  cross  section  is  obtained  directly  from 
Hansen's  rate  coefficient  by  taking  the  inverse  Laplace  transform  [68]. 
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The  zero  cross  section  for  the  recombination  of  electrons  with  an  energy 
greater  than  w  may  be  due  to  a  sudden  increase  in  the  inelastic  scatter- 
ing cross  section  due  to  vibrational  excitation  if  w = hv ,  the  threshold 
for  vibrational  excitation  by  electron  impact. 

The  two  cross  sections  are  displayed  in  Figure  5  and  the 
average  energy  of  recombining  electrons  is  shown  in  Figure  6  for  each 
cross  section.   It  should  be  noted  that  energy  transfer  due  to  inelastic 
scattering  will  be  significant  for  the  second  cross  section  at  high 
electron  energies,  which  violates  one  of  the  assumptions  made  in  deriv- 
ing the  formula  for  electron  energy  change.   Therefore,  the  effect  on 
the  electron  temperature  will  not  be  as  dramatic  as  might  be  supposed 
from  Figure  6.   It  would  be  necessary  to  add  a  term  for  inelastic 
scattering  in  Equation  (38)  in  order  to  calculate  the  variation  of  T^ 
as  a  function  of  time  for  this  cross  section. 

It  can  be  seen  from  Figure  6    that  the  recombination  heating 
effect,  due  to  the  energy  difference,  Cg-€j^^>  between  the  recombining 
electrons  and  the  average  electron  energy,  will  be  more  pronounced  for 
cross  section  a      .      The  difference  between  the  two  cross  sections  will 
be  small  at  low  temperatures  and  will  increase  as  the  temperature 
increases.   Cross  section  a  will  yield  substantially  the  same  thermal- 
ization  rates  as  the  cross  section  ct=CG   for  different  atom  temper- 
atures, but  as  the  vibrational  excitation  of  the  molecular  ions 
increases,  a  will  give  a  smaller  recombination  rate  coefficient. 

It  may  be  difficult  to  determine  the  role  of  vibrational  excita- 
tion in  the  dissociative  recombination  process.   Vibrationally  excited 
ions  may  participate  in  " superelastic"  collisions  with  electrons,  where 
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Figure  5.   Plot  of  a^  and  a^^    as  functions  of  electron  enerfry 

III  ^■' 

(arbitrary  units)  for  the  ion  ground  state  vibra- 
tional level. 
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Figure   6.       Plot   of    (f-^)/f        vs   T   for    cross    sections   of   a      and   r 
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the  vibrational  energy  is  transferred  to  kinetic  energy  of  the  electrons. 
This  would  lead  to  electron  heating,  and,  consequently,  the  observation 
of  a  smaller  recombination  rate  coefficient. 


III.   DISCUSSION  OF  PUBLISHED  EXPERIMENTAL  STUDIES 

A.   Electron  Thermalization 

Several  reports  of  electron  heating  and  long  term  electron- 
atom  temperature  disparities  in  recombining  plasmas  have  been  published. 
Various  techniques  exist  for  measuring  electron  temperatures  in  shock 
tubes  and  afterglows.   Heald  and  Wharton  [71]  have  reviewed  the  advan- 
tages and  disadvantages  of  the  several  methods  available.   It  is  suf- 
ficient to  state  that  no  one  method  gives  unambiguous  results  at  all 
electron  temperatures  and  densities  and  that  two  methods  may  give  con- 
flicting results  in  the  same  system.   Measurements  may  also  be  influ- 
enced by  changes  in  the  electron  energy  distribution,  as  noted  by 
Gregory  and  co-workers  [72] .   The  Druyvestyn  distribution 

f(Ode  =  1^1.04  6^/^/(kT)^^^J  exp  |^-0.55(6/kT)^J  d^        (47) 

has  best  been  able  to  describe  electron  energies  in  rare  gas  plasmas, 
although  the  Maxwellian  distribution  of  electron  energies  has  been 
found  to  be  satisfactory  for  plasmas  resulting  from  atmospheric  gases 
[73]. 

Many  of  the  measurements  of  electron  temperature  in  recombining 
plasmas  suggest  that  some  sort  of  electron  hefiting  is  taking  place  and 
that  the  electrons  are  not  thermal izing  as  rapidly  as  would  be  expected 
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from  the  simple  momentirr.  transfer  model  of  Petschek  and  Byron  [56]. 
Smith  and  Goodall  [3G]  -easured  the  electron  temperature  in  the  after- 
glow of  a  recombinir.j  helium-oxygen  plasma  initiated  by  a  radiofrc- 
quency  discharge.   They  noted  that  the  electron  temperature  became 
constant  at  a  temperature  about  200°K  higher  than  the  gas  temperature. 

DeLoche  et  al .  [48]  observed  elevated  temperatures  in  a  helium  after- 

17   -3 
glow  foi'  electron  densities  greater  than  10   m   .   Similar  results 

were  obtained  in  heli-o.-.  by  Fugol  et  al.  [74,75  ].   In  all  cases, 
C  »  1,  using  the  values  for  the  momentum  transfer  cross  section  of 
electrons  in  heliur.  reported  by  Crompton  et  al .  [62],  so  that  recombi- 
nation heating  should  be  negligible.   Either  some  other  mechanism  must 
be  responsible  for  The  slov-'  tliermalization  times  observed,  or  the  avail- 
able momentum  transfer  cross  sections  have  been  overestimated. 

Chen  [76]  has  r.ade  some  verj'  interesting  measurements  of  the 
electron  temperature  as  a  function  of  time  in  a  recombining  helium 
plasma.   He  observe!  that  the  electron  temperature  decreased  from 
20,000°K  to  3600°K  at  a  natural  gas  temperature  of  410°K  and  from  about 
4S,000°K  to  4200^K  at  a  neutral  gas  temperature  of  1250°K.   The  atom 

temperature  remained  constant  during  the  observation  time  (about  5  ms) , 

21   -3         ,  18   -3 
while  the  electron  density  decayed  from  4  x  10   m   to  9  x  10   m 

19   -3     ^    _18   -3 
in  the  instance  where  T  =  410°K  and  from  2  x  10   m   to  1  x  10   m 

for  the  case  where  T.  =  1250''K.   Chen  achieved  a  high  degree  of  ioniza- 

-20 

tion  and  therefore  a  lo-  value  of  C.   If  we  select  a  value  of  6.0  x  10 

r> 

m"  for  a   ,  and  substitute  Chen's  experimental  values  into  Equation  (40), 

we  obtain  C  =  0.21  for  T   =  410°K  and  C  =  29  for  T^^  =  1250' K.   This 
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would  lead  us  to  expect  slow  thermalization  in  one  case  and  more  rapid 
thermalization  in  the  latter  instance.   This  is  qualitatively  correct, 
as  is  shown  by  Figures  7a  and  7b,  which  display  a  summary  of  Chen's 
experimental  data.   However,  recombination  heating  alone  does  not 
quantitatively  explain  these  results.   In  particular,  it  does  not 
explain  a  fairly  constant  T  '^  T  ,  especially  for  T  =  1250°K. 

Electron  thermalization  has  also  been  studied  in  atmospheric 
gases.   Chen  and  Goodings  [77]  obser\^ed  the  electron  temperature  and 
density  with  electrostatic  probes  in  the  nitrogen  pinl:  afterglow  in 
a  flow  system.   They  found  that  the  electron  temperature  and  density 
increased  to  a  maximum,  then  decreased.   A  similar  time  dependence  of 
the  electron  temperature  was  found  by  Carlson  and  Rieper  [78]  in  an 
investigation  of  relaxation  phenomena  behind  shock  waves  in  nitrogen. 
Again  C  »  1,  and  recombination  heating  should  not  be  important. 

B.   Dissociative  Recombination 

One  other  explanation  for  the  slow  thermalization  of  electrons 
in  the  systems  discussed  above  is  the  presence  of  an  ionizing  source 
producing  non-thermal  electrons.   Direct  experimental  evidence  for  an 
ionizing  source  is  provided  by  measurements  of  unexpectedly  slow  rates 
of  dissociative  recombination  in  some  systems. 

Recently,  Bardsley  and  Biondi  [79]  reviewed  the  status  of 
experimental  work  on  dissociative  recombination  of  rare  gas  ions  and 
of  atmospheric  ions.   This  work  has  been  carried  out,  for  the  most 
part,  in  either  shock  tubes  or  afterglows  initiated  by  various 
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Figure  7a.   Plot  of  T^  vs  t,  as  measured  by  Chen  [76]  for  He  plasmas 


at  T   =  410°K  and  T„  =  1250°K. 
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Figure  7b.   Plot  of  a-^   as  a  function  of  time  as  measured  by  Chen  [76] 
in  a  recoinbining  He  plasma  at  T. 
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410°K  and  T  = 1250°K. 
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means.   Most  workers  agree  that  dissociative  recombination  occurs  in 
both  rare  gas  and  atmosijheric  gas  plasmas  and  that  it  has  a  large 
rate  coefficient  at  room  temperature  (300^K).   They  also  agree  that 
the  rate  coefficient  decreases  with  increasing  temperature,  but  there 
are  some  differences  of  opinion  as  to  the  actual  temperature  depen- 
dence.  Dissociative  recombination  is  known  to  yield  electronically- 
excited  atoms  in  many  instances,  and  this  fact  has  been  used  diag- 
nostically  by  Frommliold  and  Biondi  [80]  to  measure  the  rate  of  dis- 
sociative recombination  in  neon  and  argon  afterglows.   Zipf  [44]  has 
investigated  the  product  energy  distribution  in  the  recombination  of 
O   by  measuring  the  rate  constants  for  the  production  of  different 

electronic  states  of  the  oxygen  atoms.   His  data  indicate  that  0  atoms 

3      11 
are  produced  in  a  ratio  0(  P)  :  0(  D)  ;  0(  S)  :=  1.0:0.9:0.1  at  a 

temperature  of  300°K. 

The  study  of  dissociative  recombination  in  rare  gases  is 
complicated  by  the  presence  of  metastable  rare  gas  atoms.   These  may 
cause  Penning  or  associative  ionization.   The  production  of  additional 
ions  and  electrons  of  either  high  or  low  energies  could  lead  to  erro- 
neous conclusions  if  not  properly  accounted  for.   Gerardo  and  Johnson 
[81-83]  have  recently  shown  that  ion  production  is  taking  place  in 
helium  afterglow  plasmas  and  that,  contrary  to  popular  belief,  disso- 
ciative recombination  may  be  an  important  method  of  recombination  in 
helium  afterglows.   It  had  been  thought  previously  that  dissociative 
recombination  was  not  important  in  helium  because  the  m.easured  rate 
coefficient  was  so  much  smaller  than  that  observed  for  dissociative 
recombination  in  other  rare  gases. 
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Gerardo  and  Johnson  raised  the  electron  temperature  in  the 
helium  afterglow  by  heating  with  a  microwave  source.   The  electron 
heating  effectively  stopped  the  recombination  process.   It  was  observed 
that  the  concentration  of  ions  increased  dramatically,  thus  showing 
that  there  was  an  ionization  source  which  could  compete  with  recombi- 
nation processes  in  the  afterglow. 

Gerardo  and  Johnson  speculate  that  this  source  is  ionization 
due  to  collisions  of  metastable  He  atoms  in  the  afterglow.   A  recent 
paper  presenting  evidence  against  dissociative  recombination  in  He 
plasmas  is  that  by  Boulmer  et  al.  [84].   They  agree  that  the  recom- 
bination rate  coefficient  is  much  greater  than  previously  suspected 
because  of  the  presence  of  an  ionizing  source,  but  they  maintain  that 
a  dissociative  recombination  process  does  not  explain  the  experimental 
results  as  well  as  does  a  collisional-radiative  recombination  mechan- 
ism.  Obviously,  the  experimental  techniques  need  to  be  improved  still 
further  to  resolve  this  conflict. 

Besides  renewing  the  controversy  about  the  mechanism  of 
recombination  processes  in  rare  gas  plasmas,  the  measurements  of 
Gerardo  and  Johnson  point  out  the  extent  to  which  Penning  ionization 
processes  may  influence  experimental  measurements.   Experimental  mea- 
surements of  the  dissociative  recombination  rate  coefficient  of  NO 
are  therefore  suspect  when  t)ie  measurements  are  made  in  the  presence 
of  a  rare  gas  and  the  excitation  source  has  sufficient  energy  to 
produce  metastable  rare  gas  atoms. 
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The  situation  is  someuhat  simpler  for  molecular  ions  other 
than  helium.   Lai'ge  recombination  coefficients  have  been  observed 
for  Ne,  Ar ,  Kr ,  0^,  and  N  plasmas.   Table  1  summarizes  so-e  recent 
measurements  of  recombination  coefficients.   Note  the  different 
temperature  dependences  reported,  even  for  one  sj'stem.   Hagen  [85  ] 
has  also  used  a  merged  beam  apparatus  to  observe  the  electron  energy- 
dependence  of  the  cross  section  for  the  dissociative  recombination 
of  N  .   Hagen' s  results  show  an  increased  inverse  dependence  on  energy 
as  the  relative  kinetic  energy  increases.   Unfortunately,  his  exper- 
imental uncertainties  are  greatest  for  the  energies  which  are  the 
most  interesting  chemically.   Also,  his  published  results  v.ere  calcu- 
lated with  the  approximation  that 

a^  =    o^   V  (48) 

R     R 

where  a      is  an  average  cross  section  and  v  is  the  average  relative 
R 

velocity.   Entemann  and  Herschbach  [  86  ]  have  pointed  out  the  perils 
of  such  an  approximation,  especially  when  the  cross  sectior.  has  a 
strong  energy  dependence. 

Some  of  the  measured  values  of  the  dissociative  recombination 
rate  coefficient  of  NO"*"  are  presented  in  Table  2.   It  is  apparent  that 
there  are  discrepancies  in  different  measurements,  even  at  the  same 
temperature.   These  are  probably  due  to  differences  in  measurement 
techniques,  but  may  also  be  due  to  different  initial  electron  temper- 
atures or  to  different  rates  of  electron  thermal ization  duo  to  the 
factors  discussed  in  Chapter  II. 


TABLE  1 

TEMPERVrURE  DEPENDENCE  OF  THE  DISSOCIATIVE 
RECOMBINATION  RATE  COEFFICIENT 
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d  2n  a 
■"  -   d^T 

Ion 

T/K 

Reference 

Comments 

"4 

0.49 

300-4600 

40 

Afterglow 

1.5 

900-3500 

67 

Shock  tube 

0.5 

300-1500 

76 

Afterglow 

^'^ 

1.5 

670-3500 

67 

Shock  tube 

0.8-1.5 

300-1200 

76 

Afterglow 

0.67 

300-4000 

87 

Afterglow 

-; 

1.5 

800-2500 

47 

Shock  tube 

< 

0.39 

300-5000 

39 

.'Afterglow 

0.37 

700-2700 

88 

Shock  tube 

°; 

0.70 

300-1200 

39 

Afterglow 

0.56 

1200-5000 

39 

Afterglow 

1.0 

205-690 

35 

Afterglow 

0.63 

600-2500 

89 

Shock  tube 

TABLE  2 


DISSOCIATIVE  RECOMBINATION  OF  NO 
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T/K 

H„13   -3 
a  X  10   m   s 
R 

m  - 

^R 

T 

Comments 

Reference 

196 

10.0 

1. 

2 

Lyman  g? 

32 

298 

4.6 

Afterglow 

358 

3.5 

300 

3.2 

Kr  1236  mm 
Afterglow 

90 

300 

5.0  i  2.0 

Chemionization 

91 

330 

5.0 

1. 

4  ±  . 

2 

Ionospheric 
Measurements 

92,  93 

300 

1.0 

1. 

2 

Estimated 

33 

200 

7.4  i  0.7 

1. 

0 

Lyman  a 

94 

300 

4.1  ±  0. 3 

Afterglow 

450 

3.1  ±  0.2 

300 

13.0 

Lyman  a 
Afterglow 

27 

300 

5.0 

Review 

95 

2200 

0.3  ±   0.1 

1. 

4  ±    . 

2 

Ionospheric 
Measurements 

92,  93 

2900 

0.4  ±   0.2 

1. 

5 

Shock  tube 

96 

4000 

0.2 

1. 

5 

Shock  tube 

20,  21 

2000-7000 

4     -1. 
(9.0±1.5)xl0  (T/K) 

■'l. 

5 

Shock  tube 

97 
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Different  initial  electron  temperatures  arise  from  different 
ionization  sources.   It  has  been  shown  [53]  that  associative  ioniza- 
tion may  give  rise  to  low  initial  electron  temperatures.   Photoioni- 
zation  may  result  in  either  high  or  low  electron  temperatures,  depend- 
ing on  the  wavelength  of  the  photoionizing  radiation. 

Slow  thermalization  may  not  be  an  important  factor  in  any  of 
the  published  NO"*"  recombination  data  because  ex]Derimental  conditions 
generally  were  such  that  C  >  1.   Recombination  heating  could  be 
important  in  studies  of  a  higher  electron  density  plasma,  or  where 
a  high  molecular  weight  carrier  gas  is  used.   The  heating  effect 
would  be  most  evident  at  low  neutral  temperatures. 

Figure  8  shows  the  rate  coefficients  of  Table  2  plotted  as 
a  function  of  temperature.   It  is  possible  to  fit  the  data  with  a 
single  value  of  n.   However,  it  has  been  pointed  out  by  Bortner  [51] 
and  Hansen  [69]  that  it  is  not  satisfactory  to  fit  both  high  and  low 
temperature  results  with  the  same  rate  coefficient  temperature  depen- 
dence.  It  should  be  noted  that  the  rate  coefficients  at  low  temper- 
atures are  probably  too  large  due  to  the  presence  of  (NO)   ions,  as 
shown  by  Weller  and  Biondi  [  94]. 

The  high  temperature  results  in  shock  tubes  have  been  tested 
by  use  in  models  of  re-entry  flowfields.   Huber  and  co-workers  [52] 
compared  actual  re-entry  conditions  to  the  results  predicted  from  a 
theoretical  model.   They  found  good  agreement  between  their  results 
when  Dunn  and  Lordi's  [  97]  rate  coefficient  for  NO   recombination 
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Figure  8.   The  dissociative  recombination  rate  coefiicient  as  a  function  of 
temperature.   The  numbers  refer  to  the  literature  reference. 
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was  used.   More  recently,  Kiel  and  co-v.orkers  [  98]  performed  a  similar 
comparison.   Again  the  rate  coefficient  of  Dunn  and  Lordi  gave  the 
best  fit  to  experimental  data. 


IV.   THE  EFFECT  OF  AN  IONIZING  SOURCE  ON  THE  APPARENT 
RECOMBINATION  COEFFICIENT 


Experimental  results  show  that  an  ionizing  source  is  present 
in  helium  afterglows  and  suggest  that  a  similar  source  may  be  present 
in  other  recombining  plasmas.   Such  a  source  of  electrons  could  influ- 
ence measured  recombination  coefficients  and  electron  temperatures. 
If  a  source  of  electrons  is  present  in  the  recombining  system,  the 
I'ate  law,  Equation  (9)  ,  becomes 

"^"e        2 

-J—  =  -  an     +   (7  (49) 

dt       R  e    S 

where  a   is  the  rate  of  production  of  electrons  by  the  ionizing  source. 

If  the  rate  constant  is  obtained  from  the  slope  of  the  plot  of 

electron  concenti'ation  as  a  function  of  time,  then  the  apparent  rate 

constant,  a        ,    as  defined  by 
app 

<:ln  2  , 

-—  =     -  a  n  (50) 

dt      app   e 

is  related  to  the  actual  recombination  rate  constant  and  to  the  source 
rate  by 

a         =   a^   -  -^  .         (51) 

app     R     2, 
n 
e 

To  consider  how  an  ionizing  source  may  affect  the  electron  thernal- 
ization,  Eqiiations  (21)  and  (22)  must  be  modified  to 
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48 


dt    n 


r-dE, 


dt      R   e   e    s   ej 


(52) 


and 


dE      dE  ,    dE„^   dE 
tot     el     RC     s 

+  — - —  + 


dt 


dt 


dt     dt 


(53) 


The  electron  energy  change  due  to  momentum  transfer  (dE   /dt)  and 

recombination  (dE  ^/dt)  remairs  unchanged.   However,  there  is  the 
RC 

added  term  due  to  the  ionizing  source,  which  produces  electrons  with 
an  average  energy  6  . 


dE 

dt     s   s 


(54) 


Substituting  in  Equation  (52) 


dt 


_1_ 

n 


T  3/2 

e 


] 


or,  analogous  to  Equation  (38), 


(2-n) 


(55) 


(56) 


where 


a 


^s  '      A  o2 

a   n 
R  e 


(57) 


and 


6      =   I   /G 
s    s   e 


(58) 


This  suggests  that  the  effect  of  an  ionizing  source  will  have  its 
greatest  effect  on  the  observed  rate  and  electron  temperature  at 


longer  reaction  times,  when  n   is  small. 
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A  solution  to  Equation  (56)  for  C  =0.1,   /  =  55. 

s         ^s 

C  =  0.21,  n  =  1,  and  /  =  49,  is  shown  in  Figure   9  .   This  may  be 
compared  to  Chen's  [76]  experimental  data  as  shown  in  Figure  7a. 

At  high  temperatures,  collisional  ionization  and  associative 
ionization  can  be  significant  sources  of  ionization.   At  lower  temper- 
atures. Penning  ionization  and  other  processes  involving  metastable 
atoms  appear  to  be  the  most  likely  mechanisms  for  electron  production. 
Rare  gases  in  particular  have  highly  energetic  metastable  states  and 
are  commonly  used  as  a  buffer  gas  to  reduce  diffusion  of  ions  in  a 
recombining  system. 

Two  metastable  rare  gas  atoms  always  have  enough  energy  to 

produce  ionization: 

*    *         +    - 
X  +X  =X  +X  +e 

=  X"!,  +  e~  (59) 

Such  processes  are  apparently  important  in  the  He  afterglow  [81-83]. 
Since  this  requires  collision  of  two  excited  species,  the  concentra- 
tion of  metastables  must  be  high  for  this  type  of  process  to  be 
important. 

However,  Penning  or  exothermic  associative  ionization 
mechanisms 

^  -4- 

X+A=X+A+e 

=  XA"^  +  e  (60) 

could  be  important  at  lower  metastable  concentrations.   All  the 

rare  gas  metastables,  except  Xe,  have  sufficient  energy  to  ionize  NO. 
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Figure  9.   Plot  of  (^  vs  9  for  recombination  with  an  ionizing  source 
C   =0.1,  Cr.0.21,nr.l,/   =49.,/   =56. 
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Penning  processes  have  been  observed  and  found  to  have  large  cross 
sections.   For  example,  ionization  of  NO  by  Ar  has  a  rate  which 
is  about  1/5  of  the  gas  kinetic  collision  rate  [  99]. 

A  large  metastable  population  may  be  obtained  by  collisional 
excitation  in  the  ionization  source  if  high  temperatures  are  reached. 
Significant  concentrations  of  metastables  may  also  be  excited  at 
moderate  temperatures  when  an  electromagnetic  energy  source  is  used 
initially  to  produce  the  ionization  being  studied. 


V   A  SUGGESTED  RATE  CONSTANT  FOR 
DISSOCIATIVE  RECOMBINATION 


The  principle  of  microscopic  reversibility  may  be  used  to 
relate  the  cross  section  for  dissociative  recombination  to  .hat  for 
associative  ionization.   Light,  Shuler  and  Ross  [7],  Kinsey  [100], 
.evine  [101],  Marcus  [102]  and  others  have  discussed  the  application  of 
this  princ.ple  to  bimolecular  reactions  in  general  and  have  discussed 
some  of  the  advantages  and  disadvantages  of  this  approach.   A  discus- 
sion of  the  derivation  of  the  principle  of  microscopic  reversibility 

V,   +  +>,o  nature  of  particle  collisions  may  be 
from  basic  assumptions  about  the  nature  ox  pa 

found  in  Tolman's  book  on  statistical  mechanics  [103]. 

The  general  relationship  for  the  exchange  reaction 

A  +  BC  ::  AB  +  C  ^^^^ 

is  given  by 

9   -i  — >f         2   f  ~*i 
^aHA    ^    =  ^AB^cPf  ^  (62) 

Where  i  and  f  represent  the  initial  and  final  states  of  t.e  reaction. 
.   represents  the  degeneracy  of  the  states  of  species  M,  P  is  the 
relative  linear  momentum,  and  a  represents  the  state-specified  cross 

section  for  the  reaction  proceeding  xn  the  indicated  dxrection. 

It  should  be  pointed  out  that  the  cross  sections  .vhich  are 

susceptible  to  treatment  by  the  principle  of  microscopic  reversibility 
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are  the  state-specified  cross  sections  for  formation  of  a  specific 
product  state  from  a  specific  reactant  state.   Notice  also  that  degen- 
erate states  are  assumed  to  have  the  same  cross  section — this  is  an 
assumption  which  is  justified  as  long  as  there  is  no  external  field. 

The  reaction  cross  section  for  a  given  set  of  initial  conditions 
is  the  sum  of  all  cross  sections  for  formation  of  product  states  from 
the  initial  state 

CT.(E.)  =  Z  o^^^iE.)    .  (63) 

X        X  f  1 

In  most  experimental  systems  there  is  a  distribution  of  initial  states, 
so  that  the  total  reaction  cross  section  is  given  by  averaging  over  all 
initial  states 

E  a. (E.)p  . 
a(E.)  =  i      ^      ^      ^  (64) 

where  n   represents  the  distribution  function  for  the  initial  states, 
i 

The  rate  constant  is  calculated  fron  this  cross  section  by 
averaging  over  the  distribution  function  for  initial  translational 
energy,  p^ 

03 

P       1/2 
J  a(E.)E.    p^dE. 

C(T)  = .  (65) 

;  PE-i 

If  p       is   the  Max^vell-Boltzmann   distribution   function,    then   the   temper- 

"^  E 

ature-dependent  rate  constant  is  the  Laplace  transform  of  the  energy- 
dependent  reaction  cross  section 
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The  cross  section  for  the  formation  of  a  given  product  state 
is  given  bj'  averaging  the  cross  sections  of  all  initial  states  lead- 
ing to  the  final  state 

i-*f 

E  a        (E.)p. 

a^(E.)  =  ^   ^/  '  .  (66) 

i 

1 

Smith  [104]  has  shown  how  microscopic  reversibility  leads  to 
detailed  balancing  if  Maxwell-Boltzmann  distributions  are  assumed  for 
both  reactants  and  products.   Detailed  balancing  may  be  defined  as  the 
statement  that  the  forward  and  reverse  rates  of  chemical  reaction  are 
the  same  at  equilibrium  and  that  therefore  the  ratio  of  the  forward 
and  re^'erse  rate  constants  is  equal  to  the  equilibrium  constant. 

Microscopic  reversibility  may  be  applied  to  the  associative 
ionization-dissociative  recombination  reaction 

A  +  B  =  AB^  +  e~  .  (67) 

Consider  that  the  momenta 

p.  =    (2^1. E.)^''^  (68) 

1       11 

and 

Pf  =  (2^ifE^)^^^  ■  (69) 


If  E   is  the  change  in  internal  energy  during  reaction,  then 


E.  =  E    +  e'  .  (70) 
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The  reduced  masses  nre 


and 


m.  m 

^.  .  ^  ^  (71) 

1   m^  +  Mg 


m   m 

u  =   ^^   ^   «  m   .  (72) 

f   ni,„  +  m     e 
AB    e 


The  cross  section  for  associative  ionization  is  given  by 


i-^         ""e  ^f   f-i 
a^^CE.)  =  G  -^-i  a  "-(EJ  (73) 

1        U.  E.         f 
1   1 


and  that  foz'  dissociative  recombination  by 


-,  M-.  E.   .  _ 

CT    (E„)  =   -  —  —  c         (E  )  (74) 

I     G  m   E.        1 
e   f 


whore 


ff  +   g 
NO     e 
G  =  -^^^^ ~    .  (75) 

CT      Cf 

°x  ^o 

These  formulae  can  be  used  to  check  the  validity  of  proposed 
reaction  cross  sections.   For  example,  a  simple  dissociative  recombina- 
tion cross  section  of  the  form 

a^^  (Ej)  =  JL      for  E^  >  0  (7S, 

leads  to  a  corresponding  associative  ionization  cross  section 

i-f  "e  1  * 

a. 7:    (E.)  =  GK  —  —  for  E  >  E 

(1)      1  LL.   E.  1 

1    1 

=   0         for   E.  <  E*  .  (77) 

1 
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Similarly,  a  simple  associative  ionization  cross  section  of 
the  form 


* 
'(2)  ^"i'  -  "  ^'   e" 


a^rf  (E.)  =  k'(1-^)      for  E.  >  E 


_  0  for   E.  <  E  (78) 

leads  to  an  energy  independent  dissociative  recombination  cross  section 

f-i  (E  )  -  lil  ^         for   E^  >  0  .  (79) 

"^(2)  ^  f^  -  G  m  f 

^  e 

A  more  realistic  cross  section  for  associative  ionization  is 

of  the  form 

* 
i"*^  (V.    )    =   k"    (-——)  for   E,  >  E* 

i 


a:_  (E.)  =  K 

E 


(3)  --i-       V  -2 


^  0  for   E.  <  E*  .  (80) 

This  is  suggested  by  the  classical  calculations  of  Drawin  [105], 
and  leads  to  a  dissociative  recombination  cross  section  of 

f-^i  (E  )  -     ^^  _i—  for   E^  >  0.  (81) 

(•i)         I     '-'  '"e  E  +  E 

The  general  form  of  these  two  cross  sections  is  shown  in  Figures  10 
and  11.   Reduced  cross  sections 

(a^;^  E*  cm  )/(k"  p..)  =  (X-1)/(X^  (82) 

(J)        e       1 


and 


(aj:^  E*)/(k")  .  1/X+l  (83) 


* 
are  plotted  as  fimctions  of  X  =  E/E  . 
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Figure  10.   Plot  of  reduced  cross  section  for  associative  ionization. 
Defined  by  Equation  ( 82  )  . 
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Figure  11.   Plot  of  reduced  cross  section  for  dissociative  recombi- 
nation.  Defined  by  Equation  (  83  )  . 
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The  rate  constant  coiTcsponding  to  this  cross  section  for 
dissociative  recombination  may  be  obtained  by  averaging  over  initial 
rotational,  vibrational,  and  translational  energies  of  the  ion  and 
electron.   Thus, 

CO 

Oj^  =  Z   E  /   dE^  |exp  (-E^/kT)E^gy  exp  (-V^JyJ/kT) 

V^   M^   o 


X  g^   exp  f-BM^CM^+D/kTl  a^^J  (E^) 


—  m 


1/2 


dE 


{4 


/2 


f  t^f    ^^-P  (-E,/kT);  Q^^Q^ 


(84) 


where  O   and  Q   are  the  vibrational  and  rotational  partition  functions 


% 


R 


of  NO^  and  g   and  g   the  degeneracies  of  the  vibrational  and  rotational 
V       R 

levels.   The  vibrational,  rotational,  and  translational  degeneracies 
are  canceled  by  the  degeneracy  term  in  the  cross  section,  so  that  only 
the  product  of  electronic  degeneracies,  G   ,  needs  to  be  considered. 
Substituting  the  cross  section  of  Equation  (  81)  into  Equation  (  84), 
and  canceling  degeneracy  factors, 

CO 

0/  =  Z  Z  /   dE  |exp  (-E_^/kT)E^  exp  (-V^hoJ /kT) 
^   V^M^   o 


X  exp 


-BM^(M^+l)/kT 


*   — 1 
K 'V.  (G  ,m  [E.+E  ]) 
1   el  e  f 


1/2  -1 


m 


''^f  {^f''^  ®''''  (-E^/kT)j  Q^^Qj^. 


(85) 
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The  threshold  energy,  E  ,  for  anj^  specified  state  is  the 

threshold  energy  for  the  lowest  state,  E  ,  plus  the  vibrational 
energy  of  the  state,  plus  the  rotational  energy  of  the  state,  or,  in 
the  notation  that  we  have  adopted 

E  =  E   +  yjhu  +    BAr(M^+l)  (86) 

e    f       f  f 

where  B  is  the  rotational  constant.   Thus  the  cross  section  depends  on 

the  vibrational  and  rotational  energy  levels  as  well  as  the  trans- 

lational  energy  of  the  reactants. 

Replacing  the  summation  over  M  by  an  integral  and  making  the 

substitutions 

s'    =   E      +   V.haj+   E.  (87) 

e  f  I 

M„(M„+1)   ^  M^    .  (88) 

if  f 

one  obtains 

2 

<»  0°  pexp   -(BM  /kT). 

c^j^  =   K        2       /     dE^  ;     dM^   I 2 •/    ^^-P    (-Vr.^/kT) 

V  o  o  (S     +  B?vl   ) 

X  exp    (-E^/kT)E^  (89) 


where 


K      VI. 

K      = i .  (90) 

,1/2 

m 
e 


W!      eJ/2    exp    (-E^/kT)    dE^G^^y^ 


Using   the   integral    tables    [68]    to   evaluate   the   integral,    it    is 

found   that 

2 

"   exp    (-B?.L/kT) 


J     ^ _ —  clM     =   n  Jb/s'    ex-p    (sVkT)    erfc    (^,/s'/kT)  (91) 

o     (S'/B)    +    M^ 
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and   therefore 

a^  =   -^-   S     J      clE^   {e^    exp    (-E^/kT)  exp  (-VAo-VkT)  ^-^^^(S^^^T)  erfc  (./s  VkT . 

(92) 
This    expression  may  be   rearranged   to  give 

K   n     =°  rE     exp  (-E   /kT)  ex-p  (-V    huJ/kT) 

.^  =  4-v^a7BE    JciE^I^ ^ 

V        o  S 


Vs'/kT  exp  (-s'/kT)  erfc  Vs'/klj  .      (93) 


The  term 

ys'/kT  ex-p  (-s'/kT)  erfc  VsVkT 

is  nearly  constant  for  S  >  kT  and  is  equal  to  l/v'rr  for  large  values 
of  S  /kT.   Using  this  approximation,  Equation  (93)  becomes 


K  ==     .-E   exp  (-E  /kT)  exij  (-V  huj/kT)^, 


CV^  =  ^  ^/^^tTb   Z  J"  dE^  |- 


(94) 


^"f   ° 


Now  consider  a  rate  coefficient,  a   (V   ) ,    for  each  vibrational 
level  V  .   The  total  rate  coefficient  will  then  be  given  by 

^p  =   ^  N  %^^V)  ^^^-^ 

where  N   is  the  fraction  of  molecular  ions  in  vibrational  level  V  . 

^'f 
An  ex-pression  for  o  (V  )  is 

K„         =»     rE_^   exp  (-E./kT)  exp  (-V„  h'JJ  kT)-, 
a-^^(V^)  .^  .A^kTTB  J  dE^  l-i— : 1 ; ^ }    .     (96) 
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expj 


\Vhen   this    expression   is   rewritten    in   terms   of    S    ,    it   ])ecomes 


K  t»  ^exp    (-S    /kT)(S    -E^-V'^h'J^ 

c   (V^)  =~  VnkT/B  exp  (E   /kT)/  clS '    i- 

l   +\ 
e  f 


^I^   f  "  2  -        '    ^    'e  +V  hx     '  S' 


(97) 


which  can  be  solved  to  give 
c/j^(Vj)  =—  ^/rrkT/B(kT)  exp  (-V_j,  hoo/kT) 


{- 


E+V^htUs      ,E+V-h'jU.     ,E  +V  h'OJx 


where  E   represents  the  exponential  integral. 

Equation  (98  )  may  be  evaluated  for  the  recombination  of 

a  specific  molecular  ion,  if  the  information  to  substitute  into 

Equation  (90  )  for  the  K   term  is  available.   If  AB   is  NO  ,  the 

pertinent  data  are  those  tabulated  in  Table   2.   When  ex-pressions  for 

the  partition  functions  are  substituted  in  Equation  ( 90  ) ,  it  may  be 

rewritten  as 

K  =  K  (kT)"^'^  [1-exp  (-hw/kT)]  (99) 

R     o 

where 

K    =  (4k"  ■B\i.)/{m'"^    G  ,^^)  .  (100) 

o  1     e    el 

Equation  (98  )  then  becomes 

K    exp  (-V  htu/kT) 
a^{\'^)    =  -^ ^ [1-exTD  (-hu,/kT)] 


f        E  +V.  huj      /E  +V„h(jL^     ,E  +V  ha;,  . 


where  K   =  K  tt/2  ^/B 
oo     o 
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Figure  12  shows  values  of  a-  (V)/K    as  a  function  of  T  for 

R     oo 

several  vibrational  levels  of  NO  .   The  expression  for  the  total  rate 

coefficient,  assuming  a  Maxwell -Boltzmann  distribution  of  vibrational 

levels,  is  also  shown  as  cy  /K    as  a  function  T. 

R  oo 

It  will  be  noted  that  this  approach  leads  to  different  temper- 
ature dependencies  of  the  rate  coefficient  at  high  and  low  temperatures. 
Equation  (98)  also  predicts  the  temperature  dependence  which  is 
expected  for  the  recombination  coefficient  of  each  vibrational  level. 
This  is  accomplished  without  considering  the  repulsive  states  of  the 
N  +  0  system  as  O'Malley  [70]  does. 
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Figure  12.   Plot  of  reduced  rate  coefficients  from  Equation  (98). 
The  dotted  line  is  the  total  rate  coefficient  G-p^  Kqo- 
The  numbers  indicate  rate  coefficients  aj^ivy/'K^Q   for 
vibrationa]  levels  \'  =  0,1,2,3,4,5. 


VI.   A  MODIFIED  STATISTICAL  PHASE  SPACE  THEORY 

The  statistical  phase  space  theory  was  developed  in  its  present 
form  by  Keck,  Light,  and  others  [106-112]  as  a  simple  tool  for  calculat- 
ing reaction  cross  sections  and  product  energy  distributions  for  ele- 
mentary reactions.   The  theoi-y  attempts  to  minimize  assumptions  about 
the  dynamics  of  reactions — most  of  the  calculations  can  be  carried  out, 
using  tabulated  thermodynamic  and  spectroscopic  data. 

The  basic  assumption  of  the  statistical  phase  space  theory  is 
that  atoms  or  molecules  come  together  to  form  a  "strong  coupling  com- 
plex" which  then  decomposes  into  reaction  products  or  back  into  the 
original  reactants.   The  complex  will  be  formed  if  the  reactants  have 
enough  ti-anslational  energy  to  surmount  the  centrifugal  barrier  to 
complex  formation.   Products  will  be  formed  if  the  complex  has  suffi- 
cient total  energy  and  if  the  total  angular  momentum  can  be  conserved. 
The  probability  of  dissociation  of  the  complex  to  a  particular  product 
state  (or  to  a  particular  reactant  state,  for  that  matter)  is  the  same 
for  all  states  accessible  under  the  conservation  of  total  energy  and 
angular  momentum.   It  is  totally  independent  of  the  original  state  of 
the  I'eactants,  except  for  the  conserved  quantities.   Tlie  "coupling 
complex"  is  not  necessarily  a  long-lived  collision  complex,  as  has  been 
pointed  out  by  Bernstein  et  al .  [113].   An  impulsive  collision  may 
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satisfy  the  criteria  for  formation  of  a  "coupling  complex."   A  coupling 
complex  which  will  lead  to  reaction  can  be  formed  for  certain  values 
of  initial  kinetic  energy  and  orbital  angular  momentum, 

Pechukas,  Light,  and  Rankin  [108]  give  certain  coupling  criteria 
for  the  formation  and  dissociation  of  the  "strong  coupling  complex," 
although  they  caution  that  these  criteria  should  not  be  taken  too  ser- 
iously.  Essentially,  the  coupling  criteria  establish  a  maximum  initial 
orbital  angular  momentum  for  the  formation  of  the  complex  and  restrict 
the  orbital  angular  momentum  which  may  be  imparted  to  the  products. 

The  effective  potential  between  the  reactants  is 

C.    L.(L.  +  l)h 

V(r.)  =-^^      ^  \  (102) 

r  .      2u.  r . 
1       11 

and  reaction  can  only  occur  if  the  initial  energy  E.  exceeds  the  poten- 
tial energy  "hump"  which  may  be  found  by  setting 

dV 


dr 

1 


=  0  .  (103) 


The  maximum  orbital  angular  momentum  which  can  lead  to  reaction 

is,  therefore, 

2        1/3       2/3 
L.(L.  +  l)h.  =  6•^i.  C/   (E./2)     .  (104) 

As  an  example,  the  maximum  L.  for  the  K  +  HCl  system  using 

-20 

the  values  in  Table  B2  and  an  initial  energy  of  4.8  x  10    J  is 

L.      =  216  . 

1  max 

Another  coupling  criterion  which  may  be  used  for  endothermic 
reactions  can  be  derived  from  hard  sphere  reaction  theory.   Here  the 
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kinetic  energy  along  the  line  of  centers  of  the  two  reacting  particles 


must  exceed  the  activation  energy  E  , 


* 


E,     >  E  .  (105) 

I.e. 


Then,  from  Equation  (5) 


E,     =  E.(l-b^/r^)  ^  e""  (106) 

I.e.     1        o 

1/2 
where  r   is  the  hard  sphere  collision  radius.   Since  Lh  =  (20,. E.)    , 

O  11 

the  maximum  initial  orbital  angular  momentum  which  can  produce  reaction 
for  a  given  initial  kinetic  energy  is 

(2,.)V2^  (E.-E*)l/2 

L.      =  i ^ .  (107) 

1  max  n 


As  an  example,  the  maximum  L.  for  the  K  +  HCl  reaction  with 

-20 
an  initial  kinetic  energy  of  4. 8  x  10    J  and  producing  the  lowest 

vibrational  and  rotational  state  of  KCl  is 

L      =  157  .  (108) 

1  max 

This  restriction  of  the  initial  orbital  angular  momentum  has 
been  discussed  by  Levine  et  al.  [3]  as  a  general  characteristic  of  a 
reaction  proceeding  on  any  potential  energy  surface  which  has  a  large 
potential  energy  barrier  in  the  initial  reaction  channel. 

Similar  criteria  to  those  discussed  for  reactant  orbital 
angular  momentum  also  apply  to  the  product  states,  except  that  the 
initial  kinetic  energy  would  be  replaced  by  the  final  kinetic  energy 
and  the  parameters  for  the  reactants  would  be  replaced  by  the  corre- 
sponding parameters  for  the  products. 
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The  most  attractive  feature  of  tlie  phase  space  theory-  is  its 
simplicity  and  the  ease  with  which  calculations  can  be  made.   This  is 
possible  because  the  dynamics  of  the  molecular  intei^action  are  effec- 
tively ignored   and  therefore  there  is  no  need  to  consider  the  details 
of  the  potential  energy  surface. 

It  has  instead  been  assumed  that  all  product  states  have  equal 
probability.   The  theory  based  on  this  assumption,  which  shall  be 
referred  to  as  the  unmodified  phase  space  theory,  predicts  cross  sec- 
tions which  are  too  large  and  fails  to  predict  observed  inversions  of 
product  vibrational  populations.   If  the  assumption  of  equal  probabil- 
ities of  all  accessible  states  is  modified  to  assign  each  product  state 
a  probability  on  the  basis  of  some  dynamic  factor,  there  is  reason  to 
believe  that  the  theory  could  be  brought  more  into  line  with  e>q5eri- 
mental  results.   Wong  [114]  has  shown  that  it  is  possible  to  modify 
the  theory  by  introducing  "weight  factors"  for  the  various  product 
states.   He  derived  these  factors  from  infrared  chemiluminescence  data 
[115] •   This  is  not  completely  satisfactory  theory,  of  course,  because 
it  requires  the  inclusion  of  experimental  parameters  besides  the  struc- 
tural data  which  are  readily  available — namely,  vibrational  constants, 
dissociation  energies,  reduced  masses,  and  the  like.   However,  V/ong* s 
results  are  important  in  that  he  shows  that  the  statistical  phase 
space  theory  can  be  made  more  accurate  by  modifying  one  of  its  basic 
assumptions — that  all  accessible  product  states  are  equally  probable. 

The  concept  of  a  minimum  raoinentum  transfer  pi'inciple  during 
exchange  reactions  has  been  mentioned  previously  [101,116]  ,  but  has 
not  been  developed  as  part  of  a  systematic  theory.   If  we  consider  two 
particles  in  the  center  of  mass  co-oa-dinate  system,  it  is  obvious  that 
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the  two  particles  have  momenta  whicli  are  equal  in  magnitude  and  opposite 
in  direction.   If  the  two  particles  have  internal  structure  so  that 
mass  and/or  energy  may  be  transferred  from  one  to  the  other  in  an 
inelastic  or  reactive  collision,  then  the  momenta  after  the  collision 
are  not  necessarily  equal  to  the  momenta  before  the  collision.   The 
principle  of  minimum  momentum  transfer  may  be  foi'mulated  as  follows: 
"For  given  initial  conditions,  the  most  probable  product  state  is  that 
state  which  requires  the  minimum  momentum  transfer  in  the  center  of  mass 
system." 

This  assumption  is  suggested  by  theories  of  scattering  at  high 
energies.   The  spectator  stripping  model  of  chemical  reactions  assumes 
no  transfer  of  momentum  to  the  "spectator"  atom.   The  Born  approxima- 
tion, when  applied  to  reactive  and  non-reactive  scattering,  predicts 
that  the  largest  contribution  to  the  scattering  cross  section  is  due 
to  processes  involving  a  minimum  momentum  transfer.   The  Franck-Condon 
principle  assumes  that  minimum  momentum  transfer  leads  to  the  greatest 
transition  probability  between  initial  and  final  states. 

The  Franck-Condon  principle  is  discussed  in  Herzberg' s  book  on 
diatomic  molecules  [117].   The  principle  states  that  during  an  electronic 
transition  in  a  molecule  the  relative  positions  and  velocities  of  the 
nuclei  remain  constant.   In  addition  to  being  used  extensively  in  the 
interpretation  of  optical  spectra,  the  Franck-Condon  principle  can  be 
used  to  describe  such  processes  as  dissociative  recombination  [38]  and 
Penning  ionization  [118] . 
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The  Born  approximation  is  only  valid  at  high  energies  but  can 

be  extended  to  lower  energies  with  some  success,  particularly  in  regard 

to  the  angular  distribution  of  products.   Extended  discussions  of  the 

Born  approximation  and  its  applications  to  various  scattering  problems 

are  available  in  many  books  on  scattering  processes,  including  those  by 

Mott  and  Massey  [119]  and  Levine  [101].   Levine  gives  the  differential 

cross  section  f(q)  as  the  Fourier  transform  of  the  interaction  potential 

V(R) 

— » — • 

f  (q)  °^  J  e  ri— ■  V(R)dR 

— »  — » 

where  R  represents  the  reaction  coordinates.   If  V(R)  is  spherically 

symmetric,  this  can  be  rewritten  as 

CO 

sin  qr 


f(q)  cc  J*  £iiLii£  rU(r)dr 


where  r  represents  the  internuclear  distance  between  the  colliding 
species.  We  can  use  several  different  interaction  potentials.  If 
U(r)  =  S(r-r  ),  the  Dirac  delta  function,  then 

f  (q)  °^   constant  . 

2 
If  U(r)  =  A   exp  (-a  r  ),  where  A   and  a   are  constants, 

2 

f  (q)  °=  ex-p  (-q  )  . 

If  U(r)  =  A   exp  (-a  r)/r,  where  A   and  a   are  constants, 
f(q)  °^  (a^  +  q^)~^. 

If  U(r)  =  A   exp  (-a  r) ,  where  A   and  a   are  constants, 

2  2-2 

f(q)  ^  (a^  +  q  ) 
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Mott  and  Massey  [119]  show  that  the  Born  approximation  predicts 

f (q)  ^  q~^ 
for  the  inelastic  scattering  of  electons  by  atoms  and  that  this  fits 
well  the  experimental  measurements  by  Lassettre  and  coworkers  [120] 
for  the  inelastic  scattering  of  electrons  by  He  atoms.   Lassettre 
and  coworkers  have  also  studied  the  inelastic  scattering  of  electrons 

by  diatomic  molecules  [121-123],   These  cross  sections  appear  to  be 

-2 
proportional  to  q    or  exp  (-q) ,  depending  on  the  particular  transition 

studied. 

The  assumption  of  a  minimum  momentum  transfer  principle  provides 
a  means  of  introducing  additional  constraints  on  the  system  without 
destroying  the  statistical  nature  of  the  calculation  or  requiring  con- 
sideration of  the  details  of  the  dynamics  of  the  collision  complex. 
Some  implications  of  this  assumption  are  immediately  obvious.   For 
example,  the  energy  of  reaction  in  an  exothermic  reaction  will  be  more 
likely  to  go  into  the  excitation  of  vibration  and  rotation  than  into 
kinetic  energy  of  the  products,  because  the  production  of  translation- 
ally  "hot"  products  would  require  substantial  transfer  of  momentum  in 
the  center  of  mass  system.   Similarly,  an  endothermic  reaction  will  be 
enhanced  more  by  the  presence  of  reactant  internal  energy  than  by  an 
equal  amount  of  energy  present  as  kinetic  energy  of  the  reactants. 


72 


In  discussing  the  phase  space  theory,  we  shall  follow  the  nota- 
tion used  by  Truhlar  andKuppermann  [124  J,  making  some  minor  changes. 
Electronic  states  are  designated  by  U,  vibrational  states  by  V,  rota- 
tional states  by  M,  and  translational  energies  by  E.   The  subscript  i 
is  used  to  represent  the  reactants  (initial)  state  and  the  subscript  f 
is  used  to  represent  the  product  (final)  state.   The  state-specified 

cross  section  for  reactants  in  states  U. ,  V. ,  M. ,  and  E.  yielding 

111  1 

products    in   states   U^ ,    V    ,    H    ,    and  E      is    denoted  by 

a(U^,    V^,    M^,    E^;    U.,    v.,    M.,    E.)    . 
The   total    energy   E      is   given  by   the   expressions 


E^   =    E.+    (V.    +    l/2)lnj.    +    B.M.  (M.+    1)    +    E    (U.) 
Til  1111  1 


\  "    -^f"^    ^^'f   ^    1/2)!^^   +    B^M^(M^+    1)    +    E(U^) 


(109) 


where  uj  is  the  vibrational  frequency,  B  is  the  rotational  constant, 
and  E(U  )  -E(U.)  -   AE  is  the  change  in  energy  due  to  reaction. 

Pechukas,  Light  and  Rankin  [108]  derive  an  expression  for  the 
state-specified  cross  section 

a(U^,    V^,    M^,    E^;    U . ,    V.,    M . ,    E.)    = 

L. 
1   max 

rrh^  2  P    (U  VAI  E    ;U  V  M  E   ) 

2p..E.  (2M.+1)       L.=0  iiix       1111 

111  1 

(110) 

where 

m(U  V   M  E    -E  K) 

P,(U^V^M^E^;U^V.M.E.)    =  E  (2K+1)    ,^^^-^ (111) 

A(M.L.K)  -''^'^ji'^J 

where  the  symbol  A(M.L.K)  denotes  a  summation  over  all  values  of  K  such 

that    M.-L.l  ^  K  ^  M.  +  L.,  m(U.\'  M^E^ ;E_,K)  represents  the  number  of 
11  1     1      ffffT 
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ways  in  whicli  the  strong  coupling  complex  can  dissociate  into  the 
specified  final  state,  and  N(E  ,K)  represents  the  total  number  of 
states  (including  reactant  states)  into  which  a  coupling  complex  of 
total  energy  E   and  total  angular  momentum  K  can  dissociate. 

N(IL,,K)  =  S   E   E   m(U  V  M  ,E  ;  E  K). 

Uf   ^f   ^^f 

The  density  of  states  in  phase  space  can  be  seen  most  easily 

by  reference  to  a  diagram  used  by  Pechukas,  Light  and  Rankin  [108]  and 

illustrated  in  Figure  13.   Figure  13  shows  a  plot  of  L   versus  M  for 

a  fixed  total  energy,  E  ,  and  a  particular  vibrational  state,  V  ,  or 

the  products.   The  number  of  phase  space  points  accessible  for  a  given 

value  of  K  is  determined  by  the  limiting  values  M.     and  L.     where 

t  max      I  max 

M      is  the  maximum  possible  rotational  quantum  number  for  the  partic- 
f  max 

ular  vibrational  state  and  total  energy  and  L  _  is  the  maximum  final 
orbital  angular  momentum  with  which  the  products  can  be  formed,  accord- 
ing to  the  coupling  criteria  employed.   It  can  be  seen  that 

(k-M^I  =S  L^  ^  K  +  M^  (112) 

but  that 


L^  ^  L^     ,   M_  ^  M.  (113) 

f    f  max    f    f  max 


where 


E   -  V  hu  -  AE 

M^    (r.L    +  1)  = h •  (114) 

f  max   f  max  B 

Thus,  the  calculation  consists  of  adding  the  number  of  (L^,M^)  combi- 
nations that  satisfy  the  limiting  conditions  for  each  value  of  K. 
This  summation  is  proportional  to  the  crosshatched  area  in  Figure  13. 
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Figure  13.   Phase  space  density  for  a  vibrational  level, 
The  crosshatched  area  is  proportional  to 

S  m(U^,V^,M_^,E^;  E^,K). 

«f 
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Equation  (111)  can  be  modified  to  include  the  minimum  momentum 
transfer  assumption  by  letting 

ni(U^,  V^,  M^,  E^;  E^.K)  =  g^  h(q)  (115) 

for  energetically  accessible  final  states,  and  m=0  for  the  inacces- 
sible final  states.   The  function  h(q)  is  a  function  of  the  total 
momentum  transfer  q.   In  this  equation  g^  represents  the  total  degen- 
eracy of  the  product  state  (U^ ,  V^ ,  M^ ,  E^) .   As  in  the  unmodified 
theory  . 


N(i;j,,K)  =  Z  E  E  m(U^,  V^ ,  M^ ,  E^;E^,K) 


(116) 


^f  ^f  "^^f 


The  magnitude  of  the  momentum  transfer  in  the  center  of  mass 

system  is  given  by 

q'  =  pS  p'  -  2p.p^cos  0  (117) 

Where  6  is  the  scattering  angle  in  the  center  of  mass  co-ordinate 

system. 

The  function  h(q)  should  have  the  properties 

0  ^  h(q)  S  1   .  (118) 

where 

^^°^  =  '  (119) 

hC^)  =  0. 

Several  different  functions  may  be  chosen  for  h(q).   In  the 

unmodified  phase  space  theory,  h(q)  is  the  function 

h(q)  =  1     if   E^  >  0  and  |p.-p^|  <q<P.  +  Pf 

h(q)  ^  0     othcrv.-ise. 
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To  introduce  the  constraint  of  minimum  momentum  transfer,  a  dependence 
on  q  must  be  added.   A  possible  form  for  this  would  be  the  function 

h(q)  =  exp  (-q/q  )  if   E^  >  0  and  |p^-P^!  <  Q  <  ^i^^^f 

h(q)  =  0  otherwise 

where  q   is  chosen  to  be  p.  so  that  the  function  will  display  the 
dependence  on  initial  kinetic  energy  suggested  by  the  Born  approximation. 
This  function  displays  proper  behavior  as  q  -•  ^^  and  as  q  -♦  0.   It  corre- 
sponds to  the  dependence  observed  expei'imentally  by  Lassettre  and 
coworkers  for  the  inelastic  scattering  of  electrons  by  oxj'gen  [123]. 
The  experimentally  observed  inelastic  scattering  cross  section  is 
plotted  in  Figure  14  as  a  function  of  the  momentum  transfer  q. 

The  "triangular"  summation  over  possible  values  of  K  to  obtain 
the  reaction  probabilities  P   as  a  function  of  initial  orbital  angular 
momentum  L.  is  fairly  straightforward  if  the  phase  space  densities  of 
the  accessible  product  and  reactant  states  are  known.   Indeed,  if  the 
phase  space  density  is  a  smooth  function  of  K,  this  summation  may  be 
replaced  by  (2M.+1) (2L.+1) ,  so  that  Equation  (110)  now  reads 

2    ^i  max        m(U   V  ,M   E  -E  K) 

a(U  ,V  .M  .E  ;U  ,V  ,M  .E  )=^^  ^  E^    (2L..1) ,.(e^,k)        ' 

i 

The  modified  phase  space  theory  can  be  tested  by  applying  it  to 

one  of  the  few  endothermic  reactions  studied  in  molecular  beams,  the 

reaction  of  potassium  atoms  with  HCl  to  give  KCl : 

K  +  HCl  -  KCl  +  H.  (120) 

This  reaction  has  been  treated  with  the  original  phase  space  theory  by 

Truhlar  [125].   Truhlar  not  only  calculated  total  cross  sections  for 
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Figure  14.   Plot  of  cross  section  S  versus  momentum  transfer  q  for  the 
inelastic  scattering  of  electrons  by  O2,  as  measured  by 
Lassettre  and  coworkers  [123] .   All  values  are  in  atomic 
units.   The  linearity  of  this  plot  indicates  that  In  S  "^  q. 
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different  initial  kinetic  energies  and  for  production  of  various  vibra- 
tional states  but  also  calculated  the  probability  of  reaction  as  a  func- 
tion of  initial  impact  parameter  for  the  different  initial  kinetic 
energies.   These  calculations  were  then  compared  to  the  results  of 
Ross,  Greene,  and  Moursund  [126]. 

Ross,  Greene  and  Moursund  measured  the  angular  distribution  of 
elastic  scattering  and  calculated  the  reactive  scattering  by  invoking 
the  optical  theorem.   Their  model  gives  the  probability  of  reaction  as 
a  function  of  the  potential  at  the  distance  of  closest  approach. 
Truhlar  converted  this  to  probability  of  reaction  as  a  function  of 

impact  parameter,   Truhlar' s  results  for  E.  =  5.90  kj  mole    are 

1 

plotted  in  Figure  15  along  with  the  values  calculated  from  Moursund' s 
data  and  the  results  of  the  modified  phase  space  theory  calculation. 
Truhlar' s  phase  space  theory  calculations  were  repeated  for  this  par- 
ticular energy,  but  no  attempt  was  made  to  verify  all  of  his  results. 
The  modified  phase  space  theory  will  give  a  different  vibra- 
tional energy  distribution  than  the  unmodified  theory.   Table  3  shows 
the  cross  sections  for  formation  of  several  vibrational  levels  of  KCl 
for  different  initial  kinetic  energies  of  HCl  +  K.   The  percentage 
of  vibrationally  excited  molecules  is  somewhat  less  than  predicted  by 
the  unmodified  theory.   Odiorne  and  Brooks  [127]  have  also  studied  this 
reaction  in  a  molecular  beam.   Their  results  cast  doubt  on  those 
reported  by  Ross,  Greene,  and  Moursund.   Odiorne  and  Brooks  were  able 
to  differentially  detect  KCl  and  non-reactively  scattered  K;  Ross, 
Greene,  and  Moursund  were  not.   The  total  reaction  cross  sections  fovmd  by 
Odiorne  and  Brooks  were  much  smaller  than  those  seen  by  Ross,  Greene 
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Figure  15,   Plot  of  probability  of  reaction  as  a  function  of  initial 
orbital  angular  momentum  for  K  +  HCl  =  KQ  +  H  at  a 
kinetic  energy  E^  =  5.90  kj  mole"!. 


Curve  1 
Curve  2 
Circles 


Unmodified  phase  space  theory. 

Modified  phase  space  theory. 

Data  of  Ross,  Greene,  and  Moursund  [126]. 
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TABLE  3 

MODIFIED  PHASE  SPACE  THEORY  CALCULATIONS  OF 
REACTION  CROSS  SECTION 
FOR  K  +  HCl  =  KCl  +  H 

(The  endothermicity  is  assumed  to  be  3.22  kj  mole   ) 


Cross  Sect 

,.    -20 
ion/10 

2 

m 

E  /kJ  mole 

1 

v^.o 

^f  =  ^ 

V^  =  2 

V^=.3 

V^  =  4 

Total 

4.18 

3.76 

0.00 

0.00 

0.00 

0.00 

3.76 

(7.21) 

(0.00) 

(0.00) 

(0.00) 

(0.00) 

(7.21) 

5.90 

6.39 

0.02 

0.00 

0.00 

0.00 

6.41 

(10.56) 

(0.02) 

(0.00) 

(0.00) 

(0.00) 

(10.59) 

9.58 

7.86 

2.17 

0.06 

0.00 

0.00 

10.08 

(11.28) 

(2.64) 

(0.07) 

(0.00) 

(0.00) 

(13.98) 

12.05 

8.01 

3.11 

0.73 

0.00 

0.00 

11.85 

(10.83) 

(3.65) 

(0.81) 

(0.00) 

(0.00) 

(15.31) 

14.48 

7.80 

3.68 

1.37 

0.24 

0.00 

13.08 

(10.27) 

(4.20) 

(1.47) 

(0. 25) 

(0.00) 

(16.18) 

18.69 

7.51 

4.21 

2.18 

0.92 

0.23 

15.03 

(9.50) 

(4.58) 

(2.19) 

(0.92) 

(0.24) 

(17.41) 

The  values  in  parentheses  were  calculated  by  Truhlar  [125] ,  using  the 
unmodified  phase  space  theory. 
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and  Moursund.   Odiorne  and  Brooks  suggest  that  Ross,  Greene  and  Moursund 
confused  reactive  and  inelastic  scattering.   They  suggest  further  that 
the  true  endothermicity  of  the  reaction  is  not  3.22  kj  mole   but 
5.9  to  12.5  kJ  mole   .   This  latter  observation  casts  doubt  on  Truhlar' s 
calculations  as  well  as  Ross,  Greene,  and  Moursund' s  results.   Unfortu- 
nately, Odiorne  and  Brooks  did  not  measure  the  probability  of  reaction 
as  a  function  of  impact  parameter. 

Table  4  shows  several  cross  sections  calculated  using  the 
modified  phase  space  theory  for  the  K  +  HCl  reaction.   ITie  endother- 
micity was  chosen  to  be  8.3  kJ  mole   ,  since  this  is  about  the  mid- 
point of  the  range  suggested  by  Odiorne  and  Brooks.   The  initial  kinetic 
energies  ai'e  chosen  to  correspond  to  several  used  in  Table  3.   Compar- 
ison of  the  two  tables  shows  that  the  choice  of  a  value  for  the  endo- 
thermicity is  critical  in  determining  the  absolute  value  of  the  cross 
section  calculated  by  any  phase  space  theory. 

Another  method  of  measuring  the  dependence  of  reaction  proba- 
bility on  impact  parameter  has  been  suggested  by  Herschbach  [1]. 
Since  reactions  of  the  type 

M  +  HX  -^  MX  +  H  (121) 

will  necessarily  have  all  (or  most)  of  the  initial  orbital  angular 
momentum  as  product  rotation,  a  measurement  of  the  rotational  energy 
distribution  offers  a  method  of  determining  the  probability  as  a  func- 
tion of  impact  parameter.   This  could  be  done  either  in  a  molecular 
beam  apparatus  or  by  infrared  chemiluminescence  teclmiques ,  providing 
that  rapid  thermalization  of  the  product  rotational  states  can  be 
avoided.   Recent  work  [128-132]  indicates  that  this  is  now  possible. 
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TABLE  4 

MODIFIED  PHASE  SPACE  THEORY  CALCULATIONS  OF 
REACTION  CROSS  SECTIONS 
FOR  K+  HCl  =KC1  +  H 

(The  endothermicity  is  assumed  to  be  8.3  kJ  mole   ) 


Cross  Sect 

ion/: 

^-20   2 
sO     m 

E  /kJ  mole 

1 

-1 

v^  =  o 

V^  =  l 

V^  =  2 

Total 

5.90 

0.00 

0.00 

0.00 

0.00 

9.59 

2.80 

0.00 

0.00 

2.80 

12.05 

6.09 

0.23 

0.00 

6.32 

14.48 

7.75 

1.28 

0.005 

9.02 

The  exothermic  reaction 

Ba  +  CI   =  BaCl  +  CI 
has  been  studied  recently  by  Menzinger  and  Wren  [133].   They  report 
measurements  of  the  vibrational  distribution  of  the  product  BaCl  in 
the  C^TT  electronic  state.   Their  results  indicate  that  BaCl  is  formed 

with  a  statistical  vibrational  energy  distribution,  in  contrast  with 

2 
molecular  beam  measurements  which  they  quote  as  indicating  that  BaCl  (X  Z) 

is  produced  with  much  vibrational  excitation. 

The  cross  sections  for  production  of  different  vibrational 
states  of  this  reaction  can,  in  principle,  be  calculated  by  perform- 
ing the  calculations  outlined  above.   However,  the  number  of  possible 
product  states  is  now  much  greater  than  before,  and  the  long-range 
interaction  potential  constant  between  Ba  and  Cl^  and  between  BaCl  and 


83 

CI  must  be  estimated  before  the  calculation  may  be  made.   Table  B3 
shows  the  other  pertinent  data  for  the  BaCl   system. 

The  cross  sections  for  production  of  the  various  vibrational 

2 
states  of  the  ground  electronic  state  of  BaCl  (X  E)  may  be  estimated 

by  taking  the  available  phase  space  for  each  vibrational  state  propor- 
tional to  M„    X  L.     calculated  for  that  state  by  Equations  (114) 
I  max   1  max 

and  (104).   If  this  is  done  and  the  probability  of  other  products  is 

considered  to  be  small,  the  cross  section  for  the  production  of  a 

*  5/6 
given  vibrational  state  is  seen  to  be  proportional  to  (E.-E  )    , 

* 
where  E   is  the  exothermicity  of  the  reaction  plus  the  vibrational 

energy  of  the  particular  vibrational  state.   Figure  16  shows  the  rela- 
tionship between  the  phase  space  available  for  different  vibrational 

states.   Figure  17a  shows  the  relative  populations  of  the  lower  vibra- 

2 
tional  levels  of  the  X  E  state,  calculated  using  the  unmodified  and 

modified  phase  space  theory. 

The  statistical  phase  space  theory  thus  predicts  that  the 

ground  electronic  state  of  BaCl  will  be  vibrationally  excited.   The 

modified  phase  space  theory  will  predict  an  even  greater  degree  of 

vibrational  excitation  because  the  momentum  transfer  necessary  to 

form  the  lower  vibrational  states  is  greater  than  that  required  to 

produce  the  higher  vibrational  states. 

2  ^ 
The  production  of  the  excited  state,  C  TT,  of  BaCl  is  anywhere 

from  36.  kJ  mole    exothermic  to  8.4  kj  mole   endothermic,  according 
to  Menzinger  and  Wren,   A  modified  phase  space  theory  calculation  of 
the  cross  sections  for  production  of  vibrational  states  of  BaCl  pre- 
dicts a  statistical  distribution  of  product  vibrational  states  in  the 
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Figure  16.   Comparative  phase  spaces  available  to  states  1  and  2  at 
the  same  energy  and  angular  momentum. 
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2 

C  Tf  state.   The  absolute  magnitude  of  the  cross  section  for  formation 

2 

of  C  Tf  BaCl  is  small  because  of  the  large  amount  of  pliase  space  theory 

and  unmodified  phase  space  theory  predict  product  vibrational  energy 

distributions  which  are  similar  to  each  other  and  to  the  experimental 

results.   This  is  shown  in  Figure  17b. 

The  approximate  method  of  calculating  relative  vibrational 

populations  is  fairly  accurate  for  the  lower  vibrational  states  in  a 

vibrational  manifold.   However,  for  the  higher  vibrational  states,  the 

approximation  that  the  available  phase  space  is  propportional  to 

L      X  M      is  no  longer  valid,  and  a  complete  modified  phase 
f  max    f  max 

space  calculation  must  be  done  instead. 

Figures  18a  and  18b  show  the  relative  populations  of  product 
vibrational  levels  for  the  exothermic  reactions 

CI  +  DI  =  DCl  +  I 

CI  +  HI  =  HCl  +  I 
as  observed  by  Anlauf  et  al.  [115],  as  calculated  by  the  unmodified 
phase  space  theory  by  V/ong  [114]  ,  and  as  approximated  by  the  modified 
phase  space  theory.   Figure  19  shows  similar  results  for  the  exo- 
thermic reaction 

H  +  F^  HF  +  F 
which  has  been  the  subject  of  a  chemiluminescence  study  by  Pol  any i 
and  Sloan  [132].   The  initial  kinetic  energy,  E.  must  be  at  least 
100  kJ  mole    in  order  to  account  for  the  observation  of  product 
vibrational  levels  up  to  V  = 10.   These  figures  show  the  possible 
utility  of  the  modified  phase  space  theory  for  calculating  produce 
vibrational  energy  distributions.   It  can  be  seen  that,  for  exothermic 
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Figure  17a.   Relative  populations  of  the  product  vibrational  levels  for 

2 

reaction  Ba +  CI   =  BaCl  (X  E)  +  CI,  as  estimated  using  both 

the  unmodified  and  modified  phase  space  theories.   Initial 
kinetic  energy,  E.  =  20  kJ  mole 
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Figure  17b.   Relative  populations  of  the  product  vibrational  levels  for 
the  reaction  Ba  +  CI2  =^  BaCl  (C^TT)  +  CI ,  as  calculated  using 
the  unmodified  pliase  space  theory  (solid  line)  ,  as  calcu- 
lated from  the  modified  phase  space  theory  (triangles)  and 
as  measured  by  I.Ienzinger  and  Wren  [133]^  (circles) , 
Initial  kinetic  energy,  E^^  =  20  kj  mole 
icity  is  assumed  to  lie  -36  k.T  mole"  . 
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Figure  18a. 


Relative  population  of  product  vibrational  levels  foi  the 
reaction  CI +  DI =  DCl +  I .   The  triangles  are  estimated 
using  the  modified  phase  space  theory  and  an  initial 
kinetic  energy,  E,  =  5  kJ  mole-^.   The  solid  line  represents 
the  vibrationai  energy  distribution  calculated  by  Wong 
[1141  v/ith  the  unmodified  phase  space  theory.   The  circles 
are  the  experimental  results  of  -Anlauf  et  al.  [115]. 
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Figure  18b. 


Relative  population  of  product  vibrational  levels  for  the 
reaction  CI +  HI =  HCl + I.   The  triangles  are  estimated 
using  the  modified  phase  space  theory  and  an  initial 


5  kJ  mole 


-1 


The  solid  line  repre- 


kinetic  energy,  Ej^ 
sents  the  vibrational  energy  distribution  calculated  by 
Wong  [114]  with  the  unmodified  phase  space  theory.  The 
circles  are  the  experimental  results  of  Anlauf  et  al.  [115], 
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Figure  19.   Relative  population  of  product  fibrational  levels  for  the 
reaction  H+F2  =  HF+F.   The  triangles  are  estimated  using 
the  modified  phase  space  theory  and  an  initial  kinetic 
energy  E^  =  100  kJ  mole"-^.   The  circles  are  the  experimental 
results  of  Polanyi  and  Sloan  [132]. 
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reactions,  the  modified  phase  space  theory  predicts  a  product  vibrational 

energy  distribution  which  is  much  closer  to  the  experimental  results 

than  the  distribution  predicted  by  the  unmodified  phase  space  theory. 

The  associative  ionization  reaction 

N  +  0  =:  NO'^  +  e~ 

can  be  treated  with  the  modified  phase  space  theory.   It  is  reasonable 

to  assume  that  all  of  the  initial  orbital  angular  momentum  goes  into 

product  rotational  angular  momentum,  and  therefore  L  = 0.   The  hard 

sphere  coupling  criteria  can  be  used  to  calculate  L. 
^  *  1  max 

Just  above  the  reaction  threshold,  the  only  possible  reaction 
product  is  the  ground  vibrational  state  of  the  ion.   However,  inelastic 
collisions  to  form  excited  states  of  N  and  0  are  also  possible,  in  addi- 
tion to  elastic  collisions.   To  a  first  approximation,  the  cross  sec- 
tion near  threshold  is 

TTh^   2 

a  =  L     P 

2E.   1  max   L 

1 

where  P   is  a  constant  times  h(q)  and  L.      comes  from  Equation  (107). 
L  1  max 

Then, 

2   ^i"^ 

a  =  Tir   constant  X  h(q). 

o   E . 

1 

Except  for  the  term  due  to  the  momentum  transfer,  this  is  seen  to  be 
just  the  hard  sphere  collision  cross  section. 

If  sufficient  energy  is  available  to  form  one  or  more  excited 
vibrational  states  of  the  ion,  it  is  obvious  that  the  cross  section  for 
the  formation  of  the  ground  vibrational  state  will  be  greater  than  the 
cross  section  for  the  formation  of  the  first  vibrational  state,  since 
the  threshold  for  formation  of  tlie  first  vibrational  state  is  higher. 
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The  calculated  cross  sections  discussed  thus  far  should  be 
regarded  as  upper  limits  to  the  cross  sections.   The  modified  phase 
space  theory  introduces  an  angular  dependence  for  reactive  scattering 
which  was  not  present  in  the  unmodified  phase  space  theory.   This  depen- 
dence should  be  greater  for  reactions  taking  place  with  large  amounts 
of  momentum  transfer  and  less  for  reactions  taking  place  with  essen- 
tially no  momentum  transfer.   It  should  be  noted  that  angular  depen- 
dence can  also  be  introduced  into  the  phase  space  theory  by  assuming 
an  asymmetric  long  range  interaction  potential  between  products  and/or 
between  reactants.   This  would  be  a  reasonable  assumption,  but  since 
even  average  interaction  potentials  have  not  been  measured  for  many 
systems,  it  would  be  difficult  to  apply  without  resort  to  further  approx- 

mations. 

The  modified  phase  space  theory  predicts  that  a  given  amount 
of  vibrational  energy  will  be  more  effective  than  a  like  amount  of 
kinetic  energy  in  enhancing  the  cross  section  of  an  endothermic  reaction. 
If  the  kinetic  energy  is  increased,  the  momentum  transfer  in  a  reactive 
collision  will  be  increased  also,  while  if  the  vibrational  energy  of 
the  reactants  is  increased,  the  momentum  transfer  in  the  reaction  may 
decrease.   This  general  behavior  was  observed  by  Odiorne,  Brooks,  and 
Kasper  [134],  in  further  studies  of  the  K  +  HCl  reaction. 

Much  work,  both  experimental  and  theoretical,  remains  to  be 
done.  In  particular,  the  momentum  transfer  function,  h(q),  needs  to 
be  defined  rigorously  and  more  measurements  of  state-specified  cross 
sections  in  different  systems  sliould  be  made. 
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APPENDIX  A.   MOMENTUM  TRANSFER  CROSS  SECTIONS 

1.   Experimental  Measurements 

Measurements  of  the  momentum  transfer  cross  sections  of 
electrons  with  atoms  and  molecules  have  been  recently  reviewed  by 
Bederson  and  Kieffer  [63]  and  Itikawa  [64].   In  general,  the  value  of 
the  momentum  transfer  cross  section  depends  on  the  electron  energy  and 
the  identity  of  the  neutral  species.   In  determining  values  to  be  used 
to  calculate  electron  thermal ization  rates,  we  need  only  consider  the 
cross  section  at  electron  energies  near  thermal  energies.   Some  of  the 
more  recent  data  for  selected  atoms  and  molecules  is  shown  in  Table  Al 
for  energies  corresponding  to  electron  temperatures  ranging  from  250° 

to  12,500°K. 

Itikawa  [64]  expresses  confidence  in  the  reported  data  for  Ng, 
but  notes  that  there  are  uncertainties  in  the  other  values.   The  cross 
section  for  electron-argon  collisions  has  recently  been  questioned  by 
Robertson  and  Rees  [65]  as  conflicting  with  their  measurements  of  the 
longitudinal  diffusion  coefficient  of  electrons  in  argon.   The  values 
for  NO  are  not  too  accurate.   Itikawa  [64]  cites  Shkarovsky  et  al .  [57] 
as  having  made  the  only  recent  (1961)  measurements.   There  are  conflicts 
among  the  values  of  Frost  and  Phelps  [60]  and  those  of  O'Malley  [58] 
and  Chen  [59]  for  Kr.   Itikawa  quotes  values  of  Braglia  et  al.  [135] 
for  Kr  and  Xe.   The  values  for  0   are  complicated  by  a  large  cross 
section  for  electron  attachment  at  low  energies  [136]. 
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TABLE  Al 

MOMENTUM  TRANSFER  CROSS  SECTIONS  OF  ELECTRONS 
WITH  ATOMS  AND  MOLECULES 


Atom  Energies        Cross  Section 

or        References         _  _ 

Molecule  (kj   mole)       (X  10   m   ) 


He  62 


Ar  60 


0.77  to  100.      5.18   to   6.85 
200   to  500.      6.99   to   6.26 


Ne  66  2.9   to  675       .469   to   2.213 


1.0   to  30        6.0    to   0.15 
30.   to  100.      0.15   to   1.0 


Kr  64,  135       2.0   to  100.      12.    to   1. 

Xe  64,  135        2.0   to  100       80.    to   1.5 


N  61  0.1   to  100       1.  5    to   10. 

2 


O  136,  137        1.0   to  100.      3.     to   8. 

2 


1.0   to  10.       10.    to   3. 
NO  57,  64         10.   to  50.       3.     to   28. 

50.   to  100.      28.    to   10. 
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2.   Theoretical  Considei'ations 

The  energy  transfer  due  to  elastic  collisions  is  given  by 
Petschek  and  Byron  [56]  by  the  general  equation 

2 

dE  .    m   ,T      , 

— -^  =  —  (—  -  1  )  n  n   r  c   (v  )v   f  (v  )  dv  (Al) 

dt    ra.  \T     /  A  e  -J   MT   e   e   e   c    e 
A   e 

where  v   is  the  electron  velocity,  ct   (v  )  the  velocity-dependent 

momentum  transfer  cross  section,  and  f  (v  )  the  distribution  function. 

e   e 

It  would  thus  be  possible  to  assume  a  velocity  dependence  for  a        and 

•^    Ml 

use  this  expression  in  the  above  equation.   The  best  results  could 
probably  be  obtained  by  fitting  the  experimental  results  to  an  equation 
of  the  form 

"^^e^  2      3 

X     +    X^v     +    \  V   +  \   V      +    .     .     .  (A2) 


dt     b     1  e    2  e    '3  e 

where  the  \' s  are  constants. 

This  approach  leads  to  an  analytical  expression.   If  the  momentum 

transfer  cross  section  is  expressed  as  a  constant  times  a  power  of  the 

kinetic  energy: 

a,^(€  )  =  c  e^  (A3) 

MT  e     me  ^ 

and  suitable  transformations  in  the  velocity  integral  in  equation  (  Al  ) 
are  made,  then  equation  (Al  )  may  be  written  as 

2  5+r  (3+^1 

dE  .   m   .T^    ,  n  n  c    (— -)        .kT  ^  ^  2  -^ 

dt     m   Vt     /    r  ^^        2\m     J 

A        e       yrr  e 
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It  should  be  noted  that  commonly  the  "effective  collision 
frequency"  is  used  to  describe  electron-atom  momentum  transfer  results. 
This  includes  the  energy  dependence  of  the  momentum  transfer  cross 
section  implicitly,  since  the  effective  collision  frequency,  (v  ^^o), 


is  defined  by  [64] 


<v    _>   =  -^  n 


^"  3^^        ^ 


m        -,5/2 
e 


2k  T 
e 


■]■ 


03 

X    r     V      CT_(v   )    exp    (-m  v  /2kT   )    dv      .  (AS) 

>J         eMTe  ee  e  e 

o 


APPENDIX  B.   MOLECULAR  PARAMETERS 


1.   N  +  0  =  NO"^  +  e 


The  nitric  oxide  ion  is  unusually  stable  to  collisional  disso- 
ciation.  It  is  isoelectronic  to  nitrogen.   Gilmore  [138]  has  recently 
reviewed  the  available  experimental  and  theoretical  data  on  the  N-0 
system  and  has  produced  a  diagram  of  the  known  potential  curves  for 
NO  and  NO"^ .   Figure  Bl  shows  a  few  of  the  electronic  energy  states  of 

NO"^  and  their  relation  to  the  ground  state  of  NO  and  to  the  energy 

-1 
states  of  the  separated  atoms.   Note  that  over  500.  kJ  mole   separate 

the  zeroth  vibrational  levels  of  the  X  Z"  and  the  a  Z""  state.   This 

means  that  for  initial  relative  kinetic  energies  of  N +  O  less  than 

about  250.  kJ  mole"^  only  one  electronic  state  of  NO"*"  is  energetically 

accessible. 

Somewhat  different  values  for  the  energy  of  the  a  E  level  have 
been  found  by  optical  spectroscopy  and  photoelectron  techniques. 
Table  Bl  shows  the  energies  (relative  to  the  ground  state  of  NO)  of 
three  electronic  states,  obtained  by  optical  [139]  means  and  by  photo- 
electron  spectroscopy  [140], 

Table  Bl  also  shows  the  vibrational  frequencies  (in  cm   )  which 
were  obtained  from  the  photoelectron  spectra  [140]. 
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TABLE  Bl 


PROPERTIES  OF  NO 


Energy 

0) 

r 

B 

State 

(kJ  mole   ) 

(kJ  mole   ) 

-10 

(10    m) 

(kJ  mole   ) 

NO  (X^TT) 

0.0 

22.78 

1.1508 

.0204 

NO"^(X''"Z'^) 

893. 

27.75 

1.0619 

.0239 

+   3  + 
NO  (a  E  ) 

1370. 

15.61 

1.06 

.024 

NO"^  (A^m 

1766. 

18.30 

1.195 

.0189 

The  equilibrium  interatomic  radius  for  the  ground  state 
of  NO  has  been  given  variously  as  1.07  x  10    m  [6]  and 
1.0619  X  10    m  [3].   Kharitonov  [141]  has  calculated  bond  lengths 
of  1.06  X  10    m  and  1.195  X  10    m  for  the  excited  states  a  Z 


and  A  IT,  respectively. 


100 


10 


m'^ix^t!') 


N(^D)  +  0(  P) 


O^AimN 


4        1 

N(  S)  +  0(  P)   Kea!>»».J-M-»wi5P 


eV 


4       3 

N(  S)  +  0(  P) 


.g;j^a»fc»>l 


sara*^ 


.  O 


R/10    m 


2.0  2.4 


Figure  Bl.   Potential  energy  curves  lor  NO  and  NO 
1  eV  =  96.44  kj  mole"^. 
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2.   K  +  HCl  =  KCl  +  H 

There  are  three  diatomic  molecules  which  may  be  formed 
from  the  KHCl  coupling  complex:   KH,  KCl,  and  HCl.   The  KH  molecule 
is  not  energetically  possible  at  any  of  the  reactant  kinetic  energies 
which  have  been  used  in  this  paper,  but  it  is  included  in  the  table 
because  it  should  be  considered  at  higher  energies.   The  values  for 
dissociation  energy,  D^,  rotational  constant  B,  and  vibrational  fre- 
quency U)  were  obtained  from  Herzberg's  book  on  diatomic  molecules  [117]. 
The  other  values  for  HCl  and  KCl  are  taken  from  a  paper  by  Truhlar  [125], 
and  the  values  for  KH  from  a  tabulation  by  Pechukas ,  Light,  and  Rankin 
[108].   Table  B2  also  includes  the  interatomic  equilibrium  distance  r^ 
and  the  hard  sphere  collision  radius  r^.   The  constant  C^  is  the  con- 
stant from  the  long  range  intermolecular  potential. 


TABLE 

B2 

PROPERTIES 

OF 

KHCl 

K+HCl 

KH+Cl 

H+KCl 

D  /kj  mole 
e 

447. 

179. 

431. 

cu/kJ  mole 

35.76 

11.8 

3.35 

-1 

B/kJ  mole 

.1266 

.0226 

.001 

-10 
r/10     m 

1.2746 

2.244 

2.79 

-10 

r  /lO     m 
o 

3.  23 

3.07 

3.07 

„  -79   6 
C  /lO     m 
i 

373 

771 

95. 

^l/10"^^  kg 

31.08 

30.99 

1.638 
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3.   Ba  +  Clg  =  BaCl  +  CI 

The  two  electronic  states  of  BaCl  that  we  are  concerned 
with  are  the  ground  (X^E)  and  excited  (C%  .   In  Table  B3 ,  they  are 
referred  to  as  X  and  C,  respectively.   The  exothermicities ,  E^,  in 
Table  B3  were  taken  from  a  paper  by  Menzinger  and  Wren  [133].   The 
other  values  come  from  Herzberg' s  book  on  diatomic  molecules  [117], 
except  for  the  BaCl  interatomic  distance  r^,  which  was  taken  from 
Sutton's  compilation  of  interatomic  distances  [142].   The  hard  sphere 
collision  radii,  r  ,  were  estimated  to  be  about  three  times  the  radius 
of  CI,  and  the  long  range  intermolecular  potential  constant  C^  was 
estimated  to  be  equivalent  to  that  measured  for  Cs  -  Cs  interactions, 
as  tabulated  by  Kramer  and  Herschbach  [143], 

TABLE  B3 
PROPERTIES  OF  BaCl„ 


Ba  +  CI       BaCl  (X)  +  CL      BaCl  (C)  +  CI 


E   /kJ   mole"^  "   ^"  -222.  -36. 

o 

OD/kJ    mole"^  «   -^^  --^    339  3.350 

-1 


B/kJ   mole 


e 


o 


1 


0.00 

-222. 

6.754 

3.339 

,0029 

.002 

1.99 

2.82 

5.4 

5.4 

7000. 

7000. 

76.905 

48.276 

.002 


r   /10-^°   m  ^    ^^  2.82  2.82 


r   /10"^°   m  ^   ^  -^-4  5.4 


C   /lO""^^   m^  '7'^^f^  7000.  7000. 


/-LQ~27  ^c^    onq  48.9,7B  48.276 
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